Compatibility, multi-brackets and 
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Abstract 

We establish an efficient compatibility criterion for a system of gen- 
eralized complete intersection type in terms of certain multi-brackets of 
differential operators. These multi-brackets generalize the higher Jacobi- 
Mayer brackets, important in the study of evolutionary equations and the 
integrability problem. We also calculate Spencer 5-cohomology of general- 
ized complete intersections and evaluate the formal functional dimension 
of the solutions space. The results are applied to establish new integration 
methods and solve several differential-geometric problems. 

Introduction and main results 

In this paper we introduce multi-brackets of non-linear vector differential 
operators. In the case of bi-brackets they coincide with the well-known Jacobi 
bracket, which is a generalization of the classical Lagrange- Jacobi bracket im- 
portant in the theory of 1st order differential equations. These latter brackets 
become the usual commutators in the case of linear differential operators and are 
widely used in mathematical physics and non-linear analysis. We apply multi- 
bracket of differential operators to establish a criterion of formal integrability 
of systems of PDEs. 

0.1. Multi-brackets of non-linear differential operators 

Our multi-bracket {Fi, . . . , Fm+i} is defined for differential operators on sec- 
tions of the trivial m-dimensional bundle over a manifold M (notice that trivi- 
alization assumption is usually implicit for commutators or bi-brackets) and its 
value is an operator of the same kind. When Fi are linear vector differential 
operators \/i : m ■ C^^{M) C^^{M), represented as rows (V,^, . . . , V™) of 
scalar linear differential operators, the multi-bracket has the form: 

m+l 

{Vi, . . . , v,„+i} = {-If mct[vi]]$i^"' ■ V,, 
fc=i 
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where Ndet is a version of non-commutative determinant and • is the product 
(one can perceive any determinant and product for a while, but we will discuss 
various versions in the sequel). For non-linear operators the bracket is obtained 
via linearization. 

If TT is a trivial vector bundle and £{F) = {£i {F), . . . , £m{F)) is a linearization 
of the operator F, then: 

{Fi, . . .,F,n+l} = ^ X! (-l)''^a(l)(-P>(l))o- ■ ■°£a{m){Fp(m)) (-F>(m+1)) • 

We will need restriction of this bracket to the system £ of PDEs Fi — 0. Let 
ord(-Fi) = l{i) be the orders of the considered operators. Denote by diff (tt, 1) 
the algebra of all (non-linear) scalar differential operators on the sections of tt. 

Define J^Fi, ...,Fk) = {VrFi \ + \t\ < s) to be the subalgebra of the 
differential ideal generated by Fi,...,Fk in diff(7r, 1), where T>r is the total 
derivative operator with respect to a multi-index r (formula in We will 

explore the following reduced multi-bracket (further discussion in t ^3.4p : 

{Fi, . . . , F,n+i} mod Ji{i)-\ — yi{m+i)~i{Fi, ■ ■ ■ , F^+i)- 

The above equivalence class can be given by other multi-brackets, more 
convenient for calculations. In the canonical coordinates (a;',p^) on the jet- 
space J°°{tt) of a rank m (vector) bundle tt over the base M with dimM = n 
the reduced multi-bracket has the following representative: 

[Fu...,F„,+^r^— J2 «gn('^)E 

to! ^ — ' ^ — ' 

creS,„+i rt 

where 

djfi, ...Jm) 

diui, . . . ,Urn) 

is the Jacobian. In other words we have (see ^3.21 for details): 

[Fi, . . . ,F„i+i]' = {Fi, . . . mod Jl(l)-\ H(rn+l)-l{Fl, ■ ■ ■ ,Fm+i). 

When TO = 1 we obtain the Mayer bracket [F, G] of scalar differential oper- 
ators. This bracket coincides with the classical Lie-Mayer bracket for the first 
order equations and is closely related to the Jacobi bracket on higher jets. We 
apphed it in the previous works ( |KLi[ [KL2I IKL3I ) to estabhsh a compatibility 
criterion for overdetermined systems of scalar equations of a certain type. 

Namely, we considered a system of complete intersection type, i.e. given 
by r < n = dim M equations which have transversal characteristic varieties 
on regular strata. In other words, the ideal generated by the symbols of the 
operators is an algebraic complete intersection. 



f^(^(T(l)i • ■ • I -P'o-(m)) 

a(pi , ... ,p™ ) 



Tl-t \-T„ 
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0. 2. Main results 

In this paper we extend the compatibility result to the case of systems of 
PDEs with arbitrary number of unknown functions. To characterize the systems 
for which the criterion is sufficient (necessity holds always) we introduce a new 
class of systems generalizing the complete intersection for the scalar case. 

The conditions informally say the system is not too overdetermined (we will 
also discuss the opposite case) and is of general kind (transversality condition). 

Definition 1. We say a system £ C J^{t^) of r differential equations on m 
unknowns is of generalized complete intersection type if 

1. m < r < n + m; 

2. The characteristic variety has dime Char^^(£) — n + m— r — 2 at each 
point Xk (z £ {we assume dim0 = —1); 

3. The characteristic sheaf K, ower Char^^ (f ) C P^T* has fibers of dimension 
1 everywhere (see i jJ. for details of the involved objects). 

The case r = m corresponds to determined systems, where the compatibility 
conditions are void, but all the statements hold for this case as well. 

The class of systems, introduced above, is included into the systems of 
Cohen-Macaulay type, introduced in IKL2I , see also the discussion of complete 
intersection for PDEs there. Note that the number r of equations, called codi- 
mension of the system £, is defined invariantly and is calculated via the Spencer 
(5-cohomology by the formula r = codini(£) — dimH*'^{£), see IKL2I . 

Define the reduced multi-bracket due to the system {Fi, . . . , Fr) by the for- 
mula 

[Fin • ■ • J Fi^ + i]£ = {Fin ■ ■ ■ 1 Fi^+i} niod J/(ii) + ...+;(i,„^i)_i(-F'l, . . . , Fr). 

Theorem A. Consider a system of PDEs 

£ = {f,(x\..., x^, u\..., u^, = I 1 < z < rj, ord(^^,) = l{i). 

1. If the system £ is formally integrable, then the multi-bracket vanishes due to 
the system, i.e. for every collection 1 < ii < ■ ■ ■ < im+i < f 

[Fi^, . . . ,Fi^^^]s = 0. 

2. Let £ be a system of PDEs of generalized complete intersection type. Then 
the system £ is formally integrable if and only if the multi-bracket vanishes due 
to the system: 

[Fi^, . . . ,^i,„_|_i]£ = 0. 
In particular, we deduce the following compatibility criterion for scalar PDEs: 
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Corollary 1. Let £ be a system of complete intersection type, i.e. given by 
r < n differential equations on one unknown function Fi[u\ = 0, . . . , Fr[u] = 
of orders Then the system £ is formally integrable iff the Mayer 

bracket vanishes due to the system, i.e. the Jacobi bracket satisfies: 

{F„ Fj} = mod Ju+i^-i{Fi, ...,Fr), \f 1 < i< j < r. □ 

Theorem]^ was amiounced in [KL4]. The corollary was established in IKL3I 
and its particular cases for n = 2 and r — 2 appeared in |KLi| and IKL2I 
respectively. We notice however that the technique used in these papers was 
different and we required an additional assumption that at least one of the 
equations has no multiple components in the characteristic variety. We remarked 
then that this condition is superfluous, but proved the claim only for the second 
order equations. Now we can remove this technical assumption completely. 

Recall fjSl lGSilGo] ) that the obstructions to integrability belong to the sec- 
ond Spencer J-cohomology group H*'^{£) (we recall the definition in ^2.\\ . Thus 
it is important to calculate this bi-graded cohomology H*'^{£) = ®iH'^'-' {£). 

Theorem B. Let £ be a system of differential equations defined by a set of 
r = codim(£) differential operators A = (Ai, . . . , A^) : C°°{ti) — > C°°{v) (can 
be of different orders). If £ is a generalized complete intersection, then the only 
non-zero Spencer d-cohomology are given by the formula: 

{TT for j — 0, 

i/ for j = 1, 

S^-'^ir* ® A'^+^-^iy for 2 < j < r + I - m {< n). 

In the above formula we describe H*'^{£) as a usual (non-graded) vector space. 
See tjS.ll for more information about grading. 

For the case of scalar systems m = dimTr = 1 we have: H*'^{£) ~ A-'V, 
< j < J'. This corresponds to the following well-known algebraic result: 
Algebra g* of codimt;* = r is a complete intersection iff its Koszul homology 
forms the exterior algebra H,{g*) ^ A' Hi{g*), < i < r f[BH)l. 

The precise obstructions to formal integrability Wi {£) are certain curvature- 
type invariants called Weyl tensors [Li] . In |KLi[ IKL2I we calculated them for 
codim(£) = 2 complete intersections in terms of Jacobi-Mayer brackets. Now 
we can generalize this result in terms of our multi-brackets: 

Corollary 2. There is a basis ei, . . . , in H*'^{£), s — {m+i) ' '^^'^ ^ bijection 
tjj between the set of power (m + 1) subsets of {1, . . ■ ,r} and the set {1, . . . , s} 
such that the graded Weyl tensor W{£) — (BiWi{£) of the system of equations 
£ = {Fi = I 1 < « < r} with l{i) — oidFi equals 

W{£)= n 

l<ii<---<im+l<r 
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This follows directly from theorems [X] and [Bj What is more interesting is 
the precise form of the basis. We calculated it for the case of 2 scalar equa- 
tions in |KLi[ IKL2I . The result immediately generalizes to arbitrary complete 
intersections. The case to > 1 is more involved and we do not discuss it here. 

Finally we give a result on the space Ss = Sols of local/formal solutions of 
the system £ C J'''(7r) of generalized complete intersection type. As before n is 
dimension of the base M of tt and m its rank. Let r be the formal codimension 
(see §2.ip of the system £, the same number as in definition [TJ 

In classical textbooks the solutions space is characterized as follows: a gen- 
eral solution (a generic point of Solg) depends on Sp functions of p variables, 
Sp-i functions of (p — 1) variables, . . . and sq constants, where S j are C artan 
characters (introduced by E. Cartan [C]; we adapt notations from jBCG'^j ). 

Here p (called genre of £) is the maximal number, such that Sp ^ 0: only 
this character Sp has absolute meaning (citing [C]). We call the number p formal 
functional dimension of the solutions space Sg and the number d = Sp formal 
functional rank. 

The above numbers are well-defined in analytic category, i. e. when the PDEs 
and the solutions are considered analytic, see Cartan's test jBCG'^ . Cartan- 



Kahler theorem guarantees integrability. For smooth equations we need to im- 
pose additional requirements on the system to ensure that the space Solg is 
non-empty and regular (see |Ho|. El [M])- In general we take p and d to be the 
formal functional dimension and rank of the space of formal solutions. 

In abstract terms the number d equals P^^ (t) , where (t) is the Hilbert 
polynomial of the symbolic module of £ and p — deg Pg (t) + 1 (see more in §5.41 
for the detailed discussion of this subject we refer to IKL5I ). 

Theorem C. Let £ be a formally integrable system of generalized complete in- 
tersection type. Denote its orders by ki, . . . , kr and the corresponding l-th sym- 
metric polynomials by Si{ki, . . . , kr) = X]n< - <i, ^ii ' ' ' ^i, . Then the space Sols 
has formal functional dimension and rank equal respectively 

p=m + n- r-l, d= Sr-rn+l{kl, • ■ • , ^r)- 

Thus in our case p is dimension of the afhne characteristic variety and when 
r = TO + n — 1 the space Sols is a d-dimensional smooth manifold. When 
m<r<m-\-n— 1 and the system £ is analytical, the space Sols is infinite- 
dimensional and a general analytic solution depends on precisely d arbitrary 
functions of to + n — r — 1 variables. 

In smooth category the above formula for the functional rank d is important 
for formulation of well-posed boundary value problems. Note also that due to 
Cauchy-Kovalevskaya theorem the above theorem holds true in the case r = to 
of determined system of PDEs. 

0.3. Discussion and plan of the paper 

The main result (Theorem \^ provides an explicit compatibility criterion. 
To our knowledge there were only two such criteria before. One is a particular 
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case of our theorem for the first order scalar systems of PDEs - this was one 
of the motivations for the appearance of the brackets (see the historical note in 
|KLi| ) and the base for Lagrange-Charpit method (see [Gouj and !j5.2p . 

Another classical result concerns the system of linear evolution equations 
and the compatibility is expressed via commutators, being thus also a special 
case of our general result. In fact, all known integrability methods use these 
simple compatibility criteria, see 

All other methods are algorithmic, but non-explicit, and are based on the 
Cartan's prolongation-projection idea. We mention two, which apply in the 
non-linear situation. One is the Spencer theory |Go|, [S] and the Weyl tensors 
in the 2"^^ (5-cohomology groups |Li| . Another uses the differential Grobner 
basis and is being implemented into computer algebra systems now jMai IHu| . 
However neither of them give precise formulas and from computational point of 
view our criterion is more effective [Ki| . 

Theorem [B] can be specified to bi-degrees, see i l5.ll This is important, since 
it yields the place, where the system becomes involutive. In fact, we think 
that the generalized complete intersections represent the class of systems, where 
the amount of prolongations to achieve involutivity is maximal. This gives a 
possibility to reduce the estimate in the Poincare (5-lemma (see |Sw) . but this 
estimate is accepted to be too large). 

Theorem [C] gives an asymptotic estimate for the Hilbert polynomial of 
the symbolic module of the system. The dimension formula is important for 
Lagrange-Charpit method of establishing exact solutions of PDEs. 

The paper is organized as follows. In Section 1 we collect the background 
on the jet-spaces and linear differential operators and establish a machinery to 
check the formal integrability. In Section 2 we review the algebraic machinery 
and develop the commutative algebra concepts finishing with a resolvent for 
generalized complete intersections. In Section 3 we introduce multi-brackets 
and discuss their properties. Non-linear differential equations are treated geo- 
metrically (as in §l.ip and we refer the reader to |KLV| for more details. 

In Section 4 we prove Theorem |A] for linear systems and then extend the 
methods to the non-linear situation. We apply our results to construct the 
compatibility complex and non- linear Spencer cohomology. Theorems |B] and [C] 
are proved in Section 5, where we also relate our results to classical integrability 
methods and multi-Poisson geometry. 

In Section 6 we apply the compatibility criterion to solve some classical prob- 
lems in differential geometry. We discuss invariant characterization of Liouville 
metrics on surfaces and the generalized Bonnet problem. Previously the com- 
patibility criterion was applied to the plane web-geometry to solve the Blashke 
conjecture and to count Abelian relations |GLi[[GL2] - This illustrates efficiency 
of our main result. 



6 



Compatibility of PDEs via multi-brackets 

1. Jet-spaces and linear differential operators 

In this sections we collect the basic knowledge of the geometric theory of 
differential equations required for our goals. 

1.1. Systems of PDEs 

Let M be a smooth n-dimensional manifold and tt : Et^ — > M a (vector) 
bundle of rank m. Two local sections si, S2 S C{^^{tt) having tangency of order 
> k are said to have the same fc-jet at x € M and the equivalence class is called 
the fc-jet Xk = 

Thus we obtain the jet- bundle nk : J'^i'^) — > M and there arc natural pro- 
jections TTfcj : J'^(7r) -/'(tt) for / < k. We denote Xi = 'Kk,i{xk)- Any smooth 
section s e C^^{'k) induces the local section ji-s : x^ [s]^ of the bundle TTfe. 

A system of PDEs £ is represented as a collection of subsets £k C J''{tt), 
k > 0, satisfying certain conditions. The first one, regularity, is that £k with 
restricted map Trk,i is a (fiber) bundle. To formulate the second condition let us 
define for a submanifold £ C J'^(7r) its i^^ prolongation by the formula 

£W = {a;fe_,_j = [s]^+'e j''+'(7r) : jks{M) is tangent to £ at Xk with order> i}. 

Thus we form f by a collection of some given equations £k and the other £i are 
obtained via the prolongation. 

So a system of different order PDEs is the following collection of subman- 
ifolds: £i — J'(7r) until a certain order Iq, at which we add some PDEs and 
get £ig C J^°{-k), then £i = £i^~''°^ for Iq < i < li, whereupon we add new 
equations, obtain 5;^, prolong this system until jet-level I2 etc. 

Following Cartan's prolongation-projection scheme we consider Tri-^s,i{£i+s) 
and if this is a proper subset of £i, the system becomes inconsistent in the sense 
that we need to add some equations not specified in the original system. 

If we wish to exclude this we obtain: The system £ is said to be compatible 

on the level A; if ^ £k+i C and Trk+i,k '■ £k+i — > £k is surjective. In the 
regular case the last map is a bimdlc projection (submersion). 

The system £ is said to be integrable to order k if it is compatible on every 
level I < k. System £ is called formally integrable if it is integrable to order 00 
(we usually assume regularity). 

We always assume there are no functional equations in £, i.e. £q = J^iir) = 
E„. The minimal I such that £1 7^ J'(7r) is called the m,inimal order Iq of the 
system. Every number I with the property £1 7^ is called an order and 

codimension of £1 in is called its multiplicity. 

Due to Cartan-Kuranishi theorem on prolongations (in the regular case) the 
set of orders ord{£) C N is finite, i.e. there exists a maximal order Zmax starting 
from which f^"'^-' = 

Every local coordinate system {x^,u^) on the bundle tt induces coordinates 
{x^,pI) on J'^(7r) (multiindex a = (ii,...,i„) has length \a\ = ^^^lis < k), 
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where pj([s]^) = ^ ^ {x). We call them canonical coordinates. 

In a sequel we will need the operator of total differential (also denoted d): 

V : C7°°(J'=W) ^ n\M) ®c^(M) C^^iJ^+^n)). 

To define D we note that every function on J*' (tt) is a scalar differential operator 
□ : C°°(7r) C°°{M) of order k. Post-composing it with a vector field 
X e V{M) we get a differential operator □' : C°°(7r) C°°{M) of order 
/c + 1, producing the needed 1^* order differential operator Vx = ix ° '■ 
C°°(J'=W) ^C°°(J^'+H7r)). 

If we write in local coordinates X = ^X^dxi, then Vx = Y^X^Vi in 
the corresponding canonical coordinates, where the operator of total derivative 
Vi = Vg^. is given by the infinite series (when applied, only finite number of 
terms act non-trivially): 

Vi = d^.+Y^pi^,^d^^. (1) 

Similarly for X € S^V{M) we get the operator r»x : C°°(J'=(7r)) ^ C°°( J*^+'(7r)). 
For instance, if cr = (ii, . . . , i„) is a multiindex, we obtain = ■ ■ ■ ■ 

1.2. Linear differential operators 

Denote by 1 the trivial one-dimensional bundle over M . Let Ak — Difffc(l, 1) 
be the C°° (M)-module of scalar linear differential operators of order < k and 
A = DkAk be the corresponding filtered algebra, Ako Ai C Ak+i- 

Consider two linear vector bundles tt, v. Denote by Diff (tt. v) = Ufe Difffe(7r, i^) 
the filtered module of all differential operators from C°°{n) to C°°{v). We have 
the natural pairing 

Difffc(/9,j/) X Diffi(7r,p) Difffe+,(7r, i/) 

given by the composition of differential operators. 

In particular, DifF(7r, 1) is a filtered left >l-module, Diff(l,7r) is a filtered 
right .4-module and they have an .4-valued .4-linear pairing 

A e Diffi(7r, 1), V e Difrfc(l,7r) ^ (A, V) = A o V e Ak+u 

with (aA, V) = a( A, V) , ( A, Va) = ( A, V)a for a e A. 

Each differential operator A : C°°{tt) C°°{v) of order I induces an A- 
homomorphism : Diff(l,7r) Diff(l,z^) by the formula: 

Difrfe(l, 7r)9V^AoVe Difrfe+;(1, u). 

Its ( , )-dual is the .4-homomorphism (j)^ : Diff (zv, 1) — > Diff (tt, 1) given by 

DiSkii^, i)9nh->noAe Difffe+i(7r, i). 
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By the very definitions of jets with J^{'k) = C°°{'Kk) we have: 

Difffe(^,t.) = HomePC (M)(^'W,C°°M), (2) 

and differential operators A are in bijective correspondence with morphisms 
: J'(7r) V via the formula A = ^^ oji, where ji : C°°(7r) -> J^'(7r) is the 
jet-section operator. 

The prolongation : j''+'(7r) of ij:^ — ipQ is conjugated to the 

^-homomorphism (j)^ : Difffc(j/, 1) — > Difffc+i(7r, 1) via isomorphism ([2]). This 
makes a geometric interpretation of the differential operator A as the bundle 
morphism. 

Similarly, the homomorphism (j)/^ : Difffc(l,7r) Difffc+;(1, j/) is conjugated 
to the following morphism: 

Hom(j'=(l),7r) Hom(j'=+'(l),z/) 

prolongation i* 

Hom(j'(j'=(l)), J'M) Hom(j'(j'=(l)),z/), 
where i : is the natural embedding. 

1.3. Systems of differential equations as modules 

A system £ of PDEs of order I associated to A is, by definition, the subbundle 
£i = Ker(V'^) C J^n). Its prolongation is £k+i = ^ = Ker(V'^) C j'^+'(7r). 
We define the dual £* as the cokernel of the operator <j)^: 

Difffc(i., 1) ^ Difffc+^TT, 1) ^ £1^1 ^ 0. 

So £* = {£*}■ We have a natural map 7r*_^i j : £* £*_^i. Then we define the 
inductive limit £'^ — linig * . Notice that it is a filtered left yl- module. Thus we 
can consider a system as a module over differential operators (I?- module). 

The dual £a = Ker(0A) C Diff(l,7r) is a right yl- module and we have 
a pairing £^ x £^ ^ A. This pairing is however non-degenerate only for 
formally-integrable systems. This follows from the following statement: 

Proposition 1. A system £ — Ker('0'^) is formally integrable iff £* are projec- 
tive C°° (M) -modules and the maps 7t*_^_i j are injective. 

Proof. The projectivity condition is equivalent to regularity (constancy of 
rank), while invjectivity of n*^^ ^ is equivalent to surjectivity of Tr^+i^i. □ 

We can associate to the above modules their symbolic analogs. Namely, 
since ST*M ig) tt = ®S^T*M (g) tt is the graded module associated to the fil- 
trated C°°(M)-module Diff(7r, 1) = UDiffj(7r, 1), the bundle morphism (j)'^ pro- 
duces the homomorphisms (symbols) : S^^^T*M ® it ^ S^T*M ® v oi our 
differential operator A. Its dual is the graded degree / morphism 

CTA : STM ®v* STM ® tt* . 
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The value (TA,a; of ua at 2: G M is a homomorphism of S'TrM-modules. 

The ST^M-module Ma = Coker(<7A,a;) is called the symbolic module at x G 
M. Its annihilator is the called characteristic ideal /(A) and the set of its zeros 
is the characteristic variety Char(A). We will always consider projectivization 
of this conical afhne variety. 

Moreover in this paper we shall complexify the symbolic module and work 
with complex characteristics. In particular, the characteristic variety becomes 
Char^(A) C P^T*M. 

Proposition 2 [Gol[S]. Forp e T*M\{Q} letm{p) C S{T.^M) = (S;>aS'T^M 
be the maximal ideal of homogeneous polynomials vanishing at p. Then localiza- 
tion (AlA)m(p) 7^ ijff the covector p is characteristic. 

The set of the localizations (A^A)m(p) 7^ for characteristic covectorsp form 
the characteristic sheaf K, over the characteristic variety Char^(A). 

If we have several differential operators A^ e Diff(7r, of different orders 
lii 1 < * < then their sum is no longer a differential operator of pure order 
A : C°°{'k) ~^ C°^{i>), V — <S)Vi. Then <f>^ is not an ^-morphism, unless we 
put certain weights to the graded components Vi. Still we have the bundle 
morphism i^}^ and the symbol map cta : STM ® v* ^ STM (g) tt* , which 
becomes a homomorphism after a suitable weighting (in i i2.3|) . This will be 
used in the next section to pursuit the theory into the general setting of various 
orders systems. 

2. Algebra of differential equations 

Here we review the basics of symbolic theory, establish the preparatory ma- 
terial and extend it to the general case of non-linear differential equations. 

2.1. Spencer cohomology 

We consider at first the symbolic theory. Let T — T^M be the tangent 
space to the base and N = Tj;gTT^^{x), xq £ Ey^, the tangent space to the 
fiber of TT. We can identify Fk{xk) = T^AT^k^k-ii^k-i)] with S'^T* (g) N and let 
9k = g{xk) — Tx^Sk n Fk{xk) be the symbol of differential equation E. Clearly 
go = N and gi = S'^T* (g) N for i < li - the minimal order of the system. 

The symbol of the de Rham operator is called Spencer (5-operator 

S : S'^T* (E)N ^ gk-irj.* ^jsi^T* 

and it maps gk to gk-i (E" T*. In other words, if 

5^^^ = {pE S-fe+iy* (E)N\Spe gk<E)T*} 

is the first prolongation, which in the regular case equals the symbol of the 
equation then gk C g'^li- 
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Definition 2. A sequence of subspaces gk C S'^T* ^ N , k > 0, is called a 
symbolic system if gk+i C g^^K 

Thus symbols of a PDE system form a symbolic system. With every such a 
system we associate its Spencer 5-complex of order k: 

-> .9fe ^ fffe-i ® T* ^ gk-2 A^T* ^ • • • ^ gk-n ® A"T* ^ 0. 

The cohomology group at the term gi (g) hPT* is denoted by H'^'^ {£, x^), though 
we usually omit reference to the point and can also write H^'^{g). 

In terms of this cohomology / is an order of the system £ if H^~^'^{£) ^ 
and multiplicity of this order is equal to m{l) = divcig^^^/ gi ~ dimiJ'~^'^(£). 

Now the number r of equations defining £, which we called in IKL2I formal 
codimension, is the total number of orders counted with multiplicities, i.e. 

oc 

r = codim(£) = ^ dinii7'=-i'i(f ). 

k=l 

By Hilbert basis theorem this number is finite. 

Consider a symbolic system g = {gj, c S^T* ^ N\k>0} and let g* = ®gl 
be its graded dual over M (or possibly C) . Then g* is an ^T-module with the 
structure operation given by the formula 

{w ■ x,p) = {m:, S^p), w e S^T, K e 5f , p e gk+u 

where 5.^ — ■ ■ ■ 6^^. iox w — ■ ■ ■ € S^T, wj G T, and (5,^;^ — i^j. o 6 : 
9t ~^ gt-i- This g* is called the symbolic module. It coincides with the module 
Ma introduced in ^1.31 in the case of linear equations of the same order. 

This module is Noetherian and the Spencer cohomology of g dualizes to the 
Koszul homology of 5*. 

2.2. Characteristic variety and Fitting ideals 

Define the characteristic ideal by I{g) = aim(g*) C ST in the ring of poly- 
nomials R = ST and the characteristic variety as the set of non-zero covectors 
V G T* such that for every k there exists a vector w G N \ {0} with v'' G gk- 
This is a punctured conical affine variety. We projectivize its complexification 
and denote it by Char"' (5) C P^T*. When g — g{xk) is the symbol of the system 
at a point Xk G £, we also denote the characteristic variety by Char!^^(£). 

Another definition of characteristic variety is given via the homogeneous 
characteristic ideal I{g) graded by the degree: / = 

Proposition 3 [S]. Char^(g) = {p G P'^T* \ f{p^) ^ Q\J f G h^k} . 

Consider the symbolic _R- module g* . Its dimension Aivcifig* is the minimal 
number d of homogeneous elements fi, . . . , fd G R oi positive degree such that 
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the quotient g* / (/i, . . . , fd)g* is a finite-dimensional vector space. Thus due to 
equality AiTHn g* = dim[i?/ ann(5*) = R/I(g)], we can interpret 

diiRfj g* = dime Char^{g) + 1 

as dimension of the affinc characteristic variety. 

A sequence of elements fi, . . . , fs S i? is called regular if fi is not a zero 
divisor in the i?-module g* / {fi, . . . , fi^i)g* . Depth of the module g* is the 
maximal number of elements in a regular sequence. The depth and dimension 
of a module g* are related by the following inequality: 

depth 5* < dim 5* 

(we shall omit sometimes the subscript R in dimn). Now g* is called a Cohen- 
Macaulay module if depth 5* — dim^* (see |BH] for details). In such a case we 
call the system £ and the corresponding symbolic system g Cohen-Macaulay. 

For an ideal I C R and an i?-module G the length of maximal G-regular se- 
quence in I is denoted depth(/, G). Depth of the ideal / is depth/ — depth(/, R) 
(this quantity is also called the grade of the ideal /). In these terms the depth 
of the module is depth g* = depth(m, f?/ Anng*), where m = (Bi>oS^T is the 
maximal ideal in R of homogeneous polynomials with positive degree. 

We shall also use codimension of the ideal /, which is defined as codim / = 
mindimi?p, the lower bound being taken over all primes p D / in i? [E^ (in 
other sources it is called height [BHIIBV] ). 

Both the depth and the codimension are geometric quantities, i.e. they are 
defined by the conical affine locus of the ideal / over C: If Rad(/) is the radical 
of /, then depth Rad(/) = depth/ and codim Rad(/) = codim/. Moreover, 
since our ring R is polynomial, for any ideal / C /? wc have the equality 

depth / — codim /. 

For a homomorphism of free /?-modules : U V with rank U = r and 
rank]/ = m denote by Jjif) the image of the map A^U (8) A^V* R induced 
by the map A^if, where * means the functor Homfl(-, R). If we choose bases for 
U and V, i.e. identify U ~ R^ and V ~ /?'", then the map (p is represented by 
an m X r matrix A and the ideal Jj (ip) C /? is generated by all j x j minors of A. 
For the case of pure order differential operator A of i)1.3l we mean: U = ST®v* , 
V = ST(g)Tr*, ip^ (7A- 

Let G = Coker((/3). By the Fitting lemma the ideal Fitti(G) = Jm-i{v) does 
not depend on representation U ^ V ^ G ^ and is called the i-th Fitting 
invariant of G. It is known that Fitto(G) C ann(G), the two terms have the same 
radicals and the equality Fitto(G) = ann(G) holds if depth ann(G) = r — m + 1 
(see PIBV]1. 

We will be interested in the dual over R map (p* : V* U* , which is the 
symbol of the collection of differential operators determining the system £. Thus 
in this case g = Ker((^*) = Coker{(p)* , whence G = g* and ann(G) = /(<?). 

Notice that the characteristic variety of the symbolic system g can be written 
via ip* : ST*(g)Tr ST*(g)iy as Char'^(.g) = {p e P^T* \ rank(p*(p) < to}, where 
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by </5*(p) : TT ^ I' WG understand the value of ip* at the covector p. Then we 
define also the characteristic sheaf (or kernel sheaf; actually it is not a sheaf, 
but just a family of vector spaces) over Char^((?) by associating to the covector 
p the subspace ICp = Ker ip* (p) C tt. 

Remark 1. The last requirement of definition [Jl dim Kp = 1 Vp E Char*'(£), 
means that the system is similar to a system of scalar PDEs and thus be treated 
via the usual Koszul complex as in IKL2 / ■ This, however, cannot he fully for- 
malized, so that we use another approach with generalized Koszul complexes. 

2.3. Application of the Buchsbaum-Rim complex 

Let Fi be some (not necessary linear) differential operators from a bundle tt to 
a bundle vi of ordF^ = l(i), 1 < i < t. Denote fi = cr{Fi) the dual symbols (i.e. 
(ta) of these operators. These fi : Ui ^ V are S'T-homomorphisms of degree 
l{i), where = ST (g) v* , V = ST (g) ?:* . Let U = ®Ui, diuiU = J2 dim = r. 

Consider the map (p* = ®/* : V* ^ U*, which is the symbol of differential 
operator F = (Fi, . . . ,Ft). Its R-dual ^p = fi : U V is a, morphism of 
i?-modules, but it is not a graded homomorphism unless the system is of pure 
order, i.e. l{i) = k. However it becomes homogeneous of degree 1 if we consider 
the weighted grading U ~ ®Ui, with the weight l{i)~^ for the i-th summand. 

We wish to find all i?-relations between the homomorphisms fi. In other 
words, we seek to determine the 1-syzygy of the map ip -.U . It is given by 
the Fitting ideal, but we better describe the whole free resolution. 

This resolution belongs to the Buchsbaum-Eisenbud family of complexes C* 
([E]), from which we are interested in the Buchsbaum-Rim complex C^: 

Here d is the multiplication by the trace element eg V ®V* C SV ® KV* [KV* 
acts on MJ via the map A(/?*), corresponding io \ E R V ® V* : 

d{a'^ (g) fei A . . . A 6f) = ^(-l)''fc((^(fej), a) a^^^ (g) 61 A ... 6,; ... A 6*. 

The sphce map e : A"+"J7 ^ A^C/ is the action of K^Lp*n e A™L/*, where 
ri e A™y* ~ i? is a volume element (generator). 

As proved in [BR] the complex is exact iff the map (p satisfies the condition 
depth Jm {(p) > r — m + 1. 

Definition 3. Let us call an R-module G generalized complete intersection if 

codimann(G') > r — m + 1 (for a minimal free resolution U ^ V ^ G ^ Q). 

Note that the usual complete intersections G = R/I satisfy this property. 

Remark 2. By the generalized principal ideal theorem of Macaulay ('E]) we 
have: codimann(G') = codimFitto(G) < r — m + 1, so that in fact we have an 
equality above. In addition, as we shall see, the module G is Cohen- Macaulay, 
whenever it is a generalized complete intersection. 
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Proposition 4. // an R-module G is a generalized complete intersection, then 
the Buchsbaum-Rim complex is exact. 

Proof. Let ip : U V he the 1-syzygy map for r — dim U, m — dim V . 
Then we have: 

depth J„i((p) — depth Fitto(G') — depth ann(G') = codimann(G) = r — m + 1, 

where the first and third equalities are general properties of Noetherian modules, 
the last one is part of the definition and the second is a property of Fitting ideal, 
mentioned in ^ 32.21 Therefore the Buchsbaum-Rim complex is exact. □ 

Remark 3. Since the polynomial ring R is an affine domain, we have \BH\ 

dim(i?/ ann(G)) = dimi? — codimann(G') = n + m — r — 1. 

Recall that a ring P is called determinental if P = S/Qs, where 5 is a regular 
Cohen-Macaulay ring and Qs is the ideal generated by s x s minors of an m x r 
matrix A such that the codimension of Qs in S is exactly (m — s -t- l)(r — s -|- 1). 
By a theorem of Eagon and Hochster such rings are Cohen-Macaulay JBVj . Let 
us also call the ideal Qs itself determinental, if this makes no confusion. 

Theorem 5. Let the symbolic module G — g* be a generalized complete inter- 
section in the sense of definition\^ and ip : ST (g) {(Bv*) — > ST tt* be the 
corresponding R-homomorphism. Then we have: 

1. The ideal Jm{^) is determinental; 

2. Fitto(G) = lig) = ann(G); 

3. G is a generalized complete intersection in the sense of definition\^ 

Proof. Let the conditions of definition [T] be satisfied. It was shown in IKL2I 
that if the characteristic sheaf K. over Char^{g) has fibers of constant dimension 
k, then codim Cha,r^{g) < I = k{r — m + k). When A: = 1 we get I = r — m + 1 
and this is exactly the codimension of the characteristic variety of g 

codim Char''((7) = r — ni + 1, 

determined by to x m minors, or equivalently by the Fitting ideal (</?). Thus 
we see that the ideal Jmif) is determinental and codim Jm{^) = r — m + I. 

This implies that the ring R/Jmif) is Cohen-Macaulay and depth J,„((/5) = 
codim Jm((^). Since Fitto(g*) = Jm{v) and I{g) = a.im{g*) have the same 
radicals we have: 

codim /(g) = codim Jm{(p) ~ depth Jm{(p) = r — m + 1. 

Thus by the results of §2.21 we conclude that Fitto(5*) — I{g)- D 

Corollary 3. The Buchsbaum-Rim complex is a resolution of the symbolic 
module g* if the latter is a generalized complete intersection. 

Remark 4. In \KL2^ we also obtained a criterion when the ideal Fitto(g*) = 
Jra{i}}) is a topological complete intersection. 
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2.4. Non-linear differential equations 

In this section we study non-linear differential equations £. A system of 
such equations can be considered as sequence of submanifolds £k C J'^{M) 
with the property D £k+i (we assume regularity, but do not require formal 
integr ability oi £). 

Let ^ = C°°{J°°tt) be the filtered algebra of smooth functions depending on 
finite jets of tt, i.e. ^ — Ui^i with — C°°(J*7r). 

Denote — C°°{£i). The projections i^i+i^i : £i+i — + £i induce the maps 
""i+i i • T?f ~^ -Si+i, so that we can form the space = Uj^f , the points of 
which are infinite sequences (/i, /j+i, . . . ) with fi € and 7r*+i ^(/i) = fi+i- 
This is a C°°(M)-algebra. If the system £ is not formally integrable, the 
set of infinite sequences can be void, and the algebra 5"^ can be trivial. To 
detect formal integrability, we investigate the finite level jets algebras S'f via 
the following algebraic approach. 

Let £ be defined by a collection F — {Fi, . . . ,Fr) of non-linear scalar differ- 
ential operators of (possibly different) orders . . . , l{r). Each Fi determines 
a sequence of smooth maps J^{'n) — > j'^^''^*^(l) and so their collection yields a 
map J^^ii^) — > J°°(i^), where v — ®Vs with dmiVs = m(s) = {^i : = s}, 
J2 Tn{s) = r. 

Pre-composition of our differential operator F : C°°{Tr) C°°{i') with other 
non- linear differential operators gives the following exact sequence of C°°{M)- 
modules 

diff (i^, 1) ^ diff (tt, 1) ^ 5^ ^ (3) 

Note that Js{Fi^ . . . , Fr) = Im(i^)s C diffs(7r, 1) is the submodule described 
in Introduction, and 

;?f -diff,(7r,l)/Ji(Fi,...,F,). (4) 

On the level of finite jets, the map F oi ^ decreases the order appropriately, 
but is not homogenous. However we can adjust this by imposing weights to the 
vector bundles as we did in Thus we can assume for simplicity that the 
operator / has pure order k. 

It is important that the terms of ^ are modules over the algebra of scalar 
"^-differential operators ^ Diff (1,1), which are total derivative operators and 
have the following form in local coordinates [KLV] : A = ^ fa'Dcr , with /o- £ 
C°°{J°°{M)). We can identify '^#'Diff(l, 1) = Uj-^ Diff^ (1, 1) with the twisted 
tensor product of the algebras ^ and Diff(l, 1) over the action 

A: for AeDiff,(l,l). 

This '^Diff(l, 1) is a non-commutative C°°(Af)-algebra. We need a more gen- 
eral 5^-module of "^-differential operators '^Diff(7r, 1) — U*^^ Diffi(7r, 1), where 

Diff,(7r, 1) = J, (^c-(M) Diff,(^, 1). 

Remark that ■^Difr(7r, 1) is a filtered <^Difl[(l, l)-module, i.e. '^Difl[j(l,l) • 
"^Diffj(7r, 1) C '^#'Diff,+j(7r, 1). 
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Define now the filtered S^^-module Difr^(7r, 1) with "^Difff (tt, 1) = (g) 
DifFi(7r, 1). Since the module Diff (tt, 1) is projective and we can identify diff (tt, 1) 
with ^5^, we have from ([4]) the following exact sequence 

0^ J'5(F) (g)Diff^(7r, 1) <^Diff,(7r,l) ^ Difff (vr, 1) ^ 0. (5) 

Similar modules can be defined for the vector bundle v and they determine 
the 5^-module £* — Uf * by the following sequence: 

■^Difff (^., 1) ^ "^Difff+fcCvr, 1) 0, (6) 

where I : diff (tt, i/) i?<E)c°°(m) Diff (tt, i^) is the operator of universal lineariza- 
tion [KLVJ . If ^i{F). 

This sequence is not exact in the usual sense, but it becomes exact in the 
following one. The space to the left is an -module, the middle term is an 
T?f+fc-inodule. The image i'_F('^#' Difff (i/, 1)) is an -module, but we generate 
by it an 5^f^j.-submodule in the middle term. With this understanding of the 
image the term 5*^^, of ^ is an 5^f_^_j,-module and the sequence is exact. In 
other words 

=<^Difff(^,l)/(yf -ImiF). 

Sequences ([6]) are nested (i.e. their union is filtered) and so we have the 
sequence 

£U - s: ^ -> 0, (7) 

which becomes exact if we treat the image of the first arrow as the corresponding 
generated S^f-module. Thus J^g* is an S^f-module with support on Eg and its 
value at a point Xs G £s is dual to the s-symbol of the system £: 

This is a geometric definition of the symbol. Equivalently we can use the 
algebraic approach as in §1.31 Graded space g — (Bgs is dual to cokernel of the 
symbol aF of F, considered as an S'T-homomorphism ST ® v* ^ ST (g) tt*, 
which depends on the point of equation £. 

Our weight-convention apply here and hence we describe the situation on 
the level of finite jets Xi € £i for a pure order k operator F, which is the case 
represented by the following exact sequence {x = T:i{xi) G M) with the dual 
symbol map: 

^ g^{x,) S'T*M (g) TT^ "i^'^ S'~''T*M (g) v^. 

Remark 5. We interpret as the algebra of all smooth functions on our equa- 
tion £. Define £* = U£* in the same manner as , taking into account that 
the map £* — > f g+i in our approach is coupled with the change of rings (^f to 
ds+i)- 'S'o we can think of £* as of sections of the symbolic bundle g over £ with 
-coefficients. Thus (note linearization in we interpret £* as the space 
ft^{£) of differential forms on £ . 
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This remark gives us a way to treat formal integrability of the system £ as 
possibility of augmenting exact sequence ([7]) with from the left (injectivity of 
£* £s+i under the change of rings). This is reduced to the question of finding 
a left resolution of complex ^ , which will boil down onto the symbolic level as 
we shall show. 

Thus linearization of the system of generalized complete intersection type 
and methods from ^2.3\ will lead to the proof of our compatibility criterion. 

3. Multi-bracket of vector differential operators 

3.1. Non-commutative determinants 

Consider the algebra A — Diff(l,l) of linear scalar differential operators. 
It is non-commutative, so no direct generalization of the determinant function 
det : A''"y^™ — > A exists (tensor product is taken over scalars, not over A). We 
view the elements of the space A™' as rows A — (Ai, . . . , Am), which act on 
columns of functions s — {si, . . . , SmY S C°°{M)"\ 

We define non-commutative determinant Ndet : A^A"^ A via the stan- 
dard formula, where order of multiplication of matrix elements corresponds to 
the order of columns: 



In other words, we define non-commutative determinant via decomposition by 
columns, i.e. if Ci{B) is the i**^ column of _B e Ma,tmxm{A) and Mij{B) is the 
minor obtained by removing row i and column j, then we have: 



Ndet(B) = Y,{-iy-^Ci{B), Ndet(Afa(S)) = ^^(-l)"-^' Ndet(M,„(B)) a„(B), 



(it is however embarrassing to write decomposition via a mid-column) . We ob- 
viously have skew-symmetry by rows and R-linearity, but we lack ^-linearity 
and skew-symmetry by columns. Thus we can write the non-commutative de- 
terminant in the form 



Note that the symbol of the non-commutative determinant is the standard 



cr(Ndet(Vi A ... A V,„)) = det(CT(Vi) A ... A cr(V,„)), 
where the symbol of order I differential operator — (V,;i, . . . , Vi„i) G ^™ 
a(V,) =(Ti(V») = (a,(V,i),...,CTKV»m)). 




m 



Ndet(Vi A ... A V™). 



determinant 
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Beware that since the components of the operator can have smaller order, it is 
not true that = (tT(Vii), . . . ,(7{Wim))- In other words, we consider the 

grading of A"^ corresponding to increasing filtration {^"'};^q- 
Denoting U* = HomR([/, y^) we have an R-linear map 

given by the formula 



S(Vi A ... A V™+i) = Ndet 



where the last column serves as a place-holder, though the result (image of S) 
we write as as a row. 

Notice that the map S(Vi A ... A Vm+i) is a right yl-homomorphism for all 
Vi e A, i.e. ImS C Ilom'\s^\A"'+\ A) C 

Remark 6. Since our constructions are algebraic they can be generalized to 
other operator algebras, like pseudo-differential operators, Fourier operators etc. 
Then multi-brackets of the next section lead to the compatibility conditions for 
the corresponding overdetermined problems. 

One can use the theory of quasi- determinants by Gelfand et al fG^RWf to 
define other multi-brackets via similar formulas. However this requires division 
and extends the class of differential operators to non-local operators. It could be 
an exciting relation between local and global aspects of compatibility. 

3.2. Multi-brackets 

At first we define multi-brackets in the linear case. 

Let T : JV" x ■ ■ ■ x ^™ Mat(m4.i)x„(^) denote the matrix formed by 

m+l 

m -\- \ vectors- rows from y^™ . Then wc define the multi- bracket 

^m+l^m ^ A!^ 

of m + 1 vector differential operators G A"^ via the operation of the last 
section and the multiplication action x Mat(m+i)xm(-4) — > A™ on 

columns of matrices: 

{Vi, . . . , V™+i} = S(Vi A ... A V™+i) • T(Vi, . . . , V™+i). 

The z*'' component of the multi-bracket is given by 

{Vi, . . . , V™+i}, = S(Vi A ... A V™+i)(a(T(Vi, . . . , V„,+i))) . 

It is easy to check that this multi-bracket coincides with the multi-bracket de- 
fined in the introduction. This bracket is skew-symmetric by its entries and 
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is K-linear. It is not however y^-lincar and docs not commute with M-linear 
transformations of A™. 

To formulate properties of this bracket we will need later an opposite multi- 
bracket {• • • }^ : A™+^^™ A"^, which is defined by the same formula except 
that the map S is changed to S^, with the place-holder in non-commutative 
determinant being put to the first column. 

Note that 

iH h^7jx+l — 17 

where Ai is the z-th subalgebra of the filtered algebra A. 

Let now Fi g diff(m -1,1), i = l,...,m-|-l, be non-linear differential 
operators of orders ord(i^i) = which we can identify with smooth functions 
of the jet-space space, i.e. elements of C°°(J'(^)(M;]R'")) c J^(J°°(M; M")). 
Then we can define the multi-bracket {Fi, . . . , Fm+i} as an operation 

A"+ij^(J°°(M;]R")) ^ J^(J°°(A/;R")) 

via the linearization operator £ : diff (to -1,1)^ C°°(J°°(M; M")) (8)c~(a/) A"", 
see [KLV] . 

Namely, exploring the formula for the linear case, we let: 

{Fi,..., F^+i} = Eilp, A ... A • T(Fi, . . . , F„,+i). 

This multi-bracket is related to the multi-bracket of linear differential operators 
via the formula 

£{F,.....,F..^,} = {eF,,...,eF,„+A- (8) 

Similarly we can define the opposite multi-bracket {■ • • }^ for non-linear dif- 
ferential operators. 

3.3. Non-commutative "Pliicker identities" 

The multi-bracket we introduced does not satisfy the Jacobi identity of 
Nambu [Nl (or generalized Poisson) multi- bracket. Neither does it satisfy the 
axioms of SH-algebras [LS] (because the background is different: If we change 
the length of the multi-bracket, the functional space changes as well). 

However there are certain properties, these brackets do satisfy. Later in this 
section we explain that they should be viewed as a kind of generalized Jacobi 
identity. For simplicity we begin with the formulation in the linear case. 

Theorem 6. Let Vi G ,4™ be linear vector differential operators, 1 < i < m + 2, 
and let Vij denote component j of V^. Then we have the identities (where 
as usual check means absence of argument) relating the multi-bracket and the 
opposite multi-bracket for 1 < i < m: 

^(-l)'={Vi, . . . , Vfe, . . . , V„+2}IVfe = ^(-l)^Vfe,,{Vi, . . . , Vfe, . . . , V,„+2}. 



19 



Compatibility of PDEs via multi-brackets 



Proof. Indeed let T,; = Ci(T) be column i of the matrix T(Vi, . . . , Vm+2), 
the map T : A™ x • • • x ^™ Mat(,„_|_2)xm(-4) being given by the same rule 

m+2 

as T (but with one more row). 

Then the right hand side of the identity is obtained by decomposing the 
determinant 



Ndct 









Vi 






( 






T 








Vm+2 







via the first and then the last column, while the left hand side of the identity 
is the result of decomposition by the last and then the first column. But these 
operations commute. Now we unite the results by j to a row. □ 

Let <T e Sm be a permutation and : A™ — > -4™ be the corresponding linear 
transformation, = (Vi,i, . . . , V^,™) 1-^ Tc,(Vi) = (Vi,£r(i), ■ ■ ■ ,'^i,a(m))- 
This action leads to conjugated multi-brackets given by 



{Vi,...,v™+i}^ = T-i{r,(Vi), 

It's easy to see that 



,T,(V™+i)}. 



{Vi, . . . , V„,+i}J„ = {T.(Vi), . . . ,T,(V„,+i)}i = {Vi, . . . , v™+i};. 



for r = 
alone: 



1 2 
m 1 



m 
m — 1 



Thus we get an identity for the multi-bracket 



Corollary 4. For the above cyclic permutation r and arbitrary vector differen- 
tial operators V,; G A"^ it holds: 

^(-i)'={Vi, . . . , Vfc, . . . , v„+2};,Vfe - ^(-i)'=Vfe,™{Vi, . . . , Vfc, . . . , v™+2}. 

We readily generalize the above formulae for non-linear operators (with the 
help of linearization operator as in the previous section). For instance, the latter 
formula becomes: 



'^i-'^)'' (j^{Fi,...,Fk,....F„ + 2}:^^k ~ ^Fk,^{Fl, ■■■,Fk,.. . ,F,„+2} ) - 0. 



Notice that this formula for m = 1 becomes the standard Jacobi identity. 
In this scalar case our multi-bracket becomes bi-bracket and it coincides with 
the classical Jacobi bracket {F, G} of scalar (non-linear) differential operators 
F,G G diff(l,l) (in the linear case F,G Q ^ it is the commutator, for the 
non- linear case see |KLVj ) . Indeed the formula is: 

J2 (^f{G, H} - e^G.H}F) = J2 {G, H}} = 0. 

cyclic cyclic 
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Thus the multi-bracket identities could be considered as generalized Ja- 
cobi identities (but not in the sense of |APB| V We however call them non- 
commutative Pliicker identities by the following reason. Consider for simplicity 
the case m = 2. 

In this case the multi-bracket of operators Vi = (Vi_i,Vi_2), i — 1,2,3 
correspond to the composition fofi of the ^-homomorphisms 



^ ^ 



A 



3 '^0^ 



A' 



where (the determinant is Ndet) 



^1 



V2,l V2,2 
V3 1 V3,2 



Vs,! V3,2 
Vi,i Vi,2 



Vi,i Vi,2 
V2,l V2,2 



Vi.i 


Vi,2 


V2,l 


V2,2 


Va,! 


V3,2 



The above sequence is not a complex (whence the multi-bracket), but its sym- 
bolic part is a complex and is actually a resolution of the module g* correspond- 
ing to the system {Vi[u] = 0, V2M = 0, = 0}. 

To perceive the properties of the multi-brackets we consider 4 vector differ- 
ential operators G y^^, 1 < i < 4, and extend the above complex to 



O^A^ 



^4 ^ ^ j2 



t>2 is the transposed matrix and 



>1 IS 



being non-commutative 



where (^0 is a 4 x 2 matrix with rows V;, 
a skew-symmetric 4x4 matrix with entries 
determinants. 

Again the symbolic sequence is exact, but the general sequence is not a 
complex and the composition 0o</'i02 gives us the desired properties of the 
multi-bracket. Clearly the above "resolution" is built on a certain determinental 
identity, which is exactly the Pliicker identity in Grassmaninan G{2,A). 

For TO > 2 we see that our non-commutative identities model the standard 
Pliicker identities in other Grassmaninans. 



3.4. Reduced brackets and coordinates 

Let £ = {Fi[u] — 0, ■ ■ ■ , Fr[u] = 0} be an overdetermined system of PDEs, 
where Fi G diff(TO • 1,1) and u = (ui,...,u„)* € C°°(M,]R™). As in the 
introduction we denote by J7s(i?) = {^A°Fi \ ord(A)-|-ordFi < s) the submodule 
generated by Fi , . . . , F,. . 

We let £* = AT/JA^) in the linear case and £* = diff^(m • 1, 1)/Jsi£) in 
the non-linear. In this way we obtain the reduced bracket 

[/i, ■ • ■ , /m+ilf = {/i, • ■ • , fm+i} mod J/_i(£:) G 

forZ = ;(/i) + --- + ^(/m+i)- 
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Remark 7. This multi-bracket appears due to the fact, that we do not have 
a unique non- commutative determinant. If we consider determinants with the 
values in the reduced (quotient) module, as it is done in the case of Dieudonne 
determinant, we will arrive to this reduced multi-bracket. 

Every fc-th order scalar differential operator G G difffe(-l,l) induces via 
linearization a map G : £* — > £*^f^- With respect to this map Theorem [6] 
implies the reduced identities: 

Theorem 7. Let £ be an over-determined system (r > m) defined by (non- 
linear) differential operators Fi, . . . , F^. Then for any subset {ii, . . . , im+2} C 
{1, . . . , r} and any j £ [1, m] we have: 

Y,{-lfF,,^,[F,,,...,h,,...,F,^^,]e^O. □ 

Notice that we do not assume integrability and so the multi-brackets occur- 
ring in the Main Theorem are not arbitrary, but vanishing of some of them gives 
certain restrictions for the rest. 

Remark 8. In general it is not true that 

{ji, ji^^, m c 

but for formally integrable £ it is. Then we can define the bracket 

[£*,Ji^{£), ■ ■ . , Ji,„(f)]e c 

Then elements € £* such that [9,Ji^{£), . . . ,Ji„X^)]£ — with respect to this 
bracket, can be interpreted as another generalization of the classical notion of 
symmetry. 

Finally we can give a coordinate representation of the introduced multi- 
bracket. For calculational purposes it is however more convenient to work with 
the following multi-bracket: 



l<i<m ^ 



m 

3 



where Fi G diffi(i)(TO • 1,1). For rn = 1 this gives Mayer brackets instead of 
Jacobi brackets |KLi| . The following statement is straightforward: 

Proposition 8. Restrictions of the two multi-brackets to the system £ coincide: 

[Fi, . . .,F,n+i]£ = [Fi,.. . mod Ji-i{£), 

where I = J2i=i^ K*) '^^ before. □ 
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4. Compatibility criterion 

In this section we prove our compatibility criterion. Its particular cases are 
theorems from |KLi[ |KL2[ IKL3I , where we used geometric theory of PDEs and 
the obstructions to compatibility were identified with certain curvatures (Weyl 
tensors). Here we propose an approach based on the construction of symbolic 
compatibility complex, which uses the dual algebraic approach. 

4.1. Syzygies for modules of linear differential operators 

Consider the filtered ^-module Diff(7r, 1) and let 

Diff (t/, 1) ^ Diff (tt, 1)^£* ^0 

be a representation of the dual to a system of linear equations £. Set I = 
Im((/)^). In other words, we let Tk C Difffe(7r, 1) denote the sequence of sub- 
modules Jk{£) = Jk{Fi, . . . , Fr) as in Introduction, Fi G Diff;(j)(7r, 1). Notice 
that these Ik define our equations: £k = {xk ■ h{xk) = OV/i e Ik} C J'^(7r). 

Proposition 9. The formal integrability of the system £ is equivalent to the 
requirement that Tk+i H Difffc(7r, 1) = 2k for all k. 

Remark 9. The last condition means it is not possible to get new relations of 
order k inX from the relations of order k + 1 via linear combinations over Ak+i ■ 

Proof. In fact, surjectivity of the projections i^k+i.k '■ £k+i ^ £k (formal 
integrability) is equivalent to injectivity of the dual maps T^k+i k • ^k ~^ ^k+i- 
The claim follows from the natural isomorphism £^ ~ Difffc(7r, 1)/Ifc. □ 

Now if T is generated by linear differential operators Fi, . . . , Fr, every ele- 
ment in Tk+i is represented in the form ^ AiFi, where Ai — o.IT>t are scalar 
differential operators of ord(Ai) < fc -I- 1 — ord(i^i). The condition of proposition 
[9] is equivalent to the following relation on the fc-symbols: 

A,F,) = J2 <^{A^)a{F,) = G S'^+^T ® N*. (9) 

To describe all such relations in the submodule of symbolic relations / C ST^N* 
we use the syzygy approach of ^'2.3\ 

We also denote V* = ST* (g) tt and U* = ST* (g) 1^, where the bundle v = ®Vi 
is, in general, graded by the degrees of operators A^. We have symbols fi = 
<j{Fi) (previously denoted a^') of our differential operators Fi, i — 1, . . . ,r. The 
map ip — (f* : V* ^ U* , which can be represented in bases as ip : i?™ , 
(ui,...,?/™) (/i,...,/'0, has the kernel g C 5r* (g iV. 

If n = n(.g) C ST<^N* is the annihilator submodule, then ST<^N*/U{g) is 
the symbolic module g* — Coker((y5). Consider R-dual to the Buchsbaum-Rim 
complex from (as the Spencer complex is M-dual to the Koszul complex): 

O^g^V* ^U* ^ A^+^C/* ^ F^® A"+2C/* ^ . . . 5"-™-iv-x ® A'-[/* ^ 0, 



23 



Compatibility of PDEs via multi-brackets 



where 5 ^ d* , lo ^ e* and = {V*)* . Choosing a basis ei, . . . , e„i of 1^* ~ 
R™- we can describe informally cj(^) = ^ A ■0(ei) A ... A ip{em)- However as 
the symbolic differential operator (p decreases degrees, we need to change it to 
. . . h,nO — ^ ^ ipihiei) A ... A ip{hmem) for certain elements hi ^ R of 
sufficiently high degrees. 

Thus all the relations between the symbolic differential operators a{Ai) are 
given by the explicit formula u o ijj = from the above complex, which we 
eventually call generalized Spencer 5-complex. 

4.2. Proof of the main theorem for hnear systems 

At first let us consider two partial cases. 

1. Scalar equations. In this case m = 1 and the condition of definition [1] says 
that g* — ST/I{g) is the usual complete intersection. So /(g) = (/i, . . . 
where {fi} form a regular sequence of length r < n. Then the Koszul complex 

-> / (g) A'R'' ^ > / «) A^R"^ ^/«)i?''->/^0 

is exact. In particular, vanishing of the 1st homology yields: 

r r 

i=i j=i 

and so the 1-syzygy module is generated by the relations fifj — fjfi ~ 0. Thus 
we need to check the condition of proposition[9]only for combinations FiFj—FjFi 
(multiplication in algebra A is non-commutative), see Q. 

As a consequence we obtain that the commutator [Fi, Fj], which is an oper- 
ator of order kj = + — 1 should belong to the space Ji-. generated by 
Fi, . . . , Fr and their total derivatives up to the order lij . This is the compati- 
bility condition for the system £, exactly as theorem lAl states. 

2. Systems on two-dimensional manifolds. In this case n — 2 and condition 
1 of definition [T] gives r = m + 1. Now instead of Koszul complex we use the 
following approximation to a resolution 

where / = aJm{''P) C i? for a non-zero divisor a. Let fi : i?'" ^ i? be symbols of 
the defining equations for the system £ and fi{ej) = fij their values on a basis. 
Denote by A{f) = \\fij\\ the {m + 1) x m matrix of the operator ip. Then the 

map r is given by the formula r(^) = J2i=i^i~^y~^^i'^^^ ^iif)' '^i^h Ai{f) is 
obtained from A{f) by deleting the i-th. row (Laplace decomposition). 

The Hilbert-Burch theorem states that the above complex is exact whenever 
depth Jm{i^) > 2. By theorem [5] this follows from the conditions of Definition [T] 
(in particular Char"' (17) = 0) and also we see that we have equality. Thus the 
only generator of 1-syzygy is the relation: 

m+l 

^ (-1)^-1 detA,(/)/, = 0. 

i=l 
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Let Ai{f) be some differential operators with the symbols Ai{f). As a conse- 
quence of the 1-syzygy description we obtain that the expression 

m+l 

^ (-1)^-1 detA(/)F„ (10) 

i=l 

which is an operator of order I = + ■ ■ ■ + l{r) ~ 1, should belong to the 
space J^i generated by Fi, . . . , and their total derivatives up to the order I. 
Noticing that the expression in (|10p is the multi-bracket {Fi, . . . ,Fr}, we get 
the compatibility condition from theorem [XI 

Now, as we have clarified the simple situations, we will study 

3. The general linear case. 

Here we should use the generalized Spencer complex described in §4.11 (R- 
dualization of the complex from ^2.3^ . If the assumptions of definition [1] are 
satisfied, we conclude that the relations generating the 1-syzygy are of the type 
(fTO)) . namely 

m+l 

Y,{-iy-' dct A,,^{f)f,^, 

i=i 

where ^^(/) is obtained from the r x m matrix A{f) of the map ip by deleting 
rows with the numbers <; C {1, . . . , r} and g is the complement to the subset 
. . . , im+i}, so that the resulting matrix Ag^if) is square of size m x m. 
It's clear that changing Ag^i^{f) i-^ Ag.ijif) and fi. i-^ Fi. in the above 
expressions we obtain the multi-brackets {Fi^ , . . . , Fi^^^ } (or [Fi^ , . . . , Fi^^^] 
depending on the manner we extend the symbol Ag,i. (/) to a differential opera- 
tor Ag.i - (/)) for various (ordered) subsets {ii, . . . , im+i\ C {1, . . . , r}, which by 

proposition [9] should belong to the subspace i7;(ii)H hK«m+i)-i ^ formally 

integrable. The claim is proved. 

4.3. The proof for non-linear operators and generalization 

Now consider the general case of non- linear systems £. The main theorem 
in this case can be proved as follows: We linearize the system. The new linear 
system is of generalized complete intersection type as well. Compatibility of 
non-linear system £ is equivalent to compatibility for each linearization on a 
jet-solution. The multi-brackets are also preserved (this is well-studied in the 
case of usual brackets, see |KLVj for a relation between universal linearization 
operators, Jacobi brackets and evolutionary differentiations), see ([5]). 

We however will not develop these vague ideas, but show instead how to 
modify the proof from ^14. 21 Let the system £ be given by a (matrix) differential 
operator A = (i^i, . . . ,F^) : C°°(7r) C°°{v). The bundle v can be locally 
trivialized (R"" above) or more generally we can split it ly — (Di^s according to 
different orders of the differential operators Fi . Since we showed in H2.'S\ how to 
reduce this case to the pure order by imposing weights on Vg, we just assume 
the operator is of pure order k (can be set 1) A : J*(7r) J^^^{v). Its symbol 
we denote by Lp. 
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Let t = i — k{m + 1). Denote by AJ^"*"^ the t-th graded component of 
the (m + l)-st exterior degree of the corresponding module in the context of 
symbolic systems and the t-th filtered submodule for the case of differential op- 
erators (with the above weight-convention) . Let us also use the short notations 
Diffj(7r) = Diff ((tt, 1), Diffj(i/) = Diffj(!/, 1) [only for large diagrams]. 

The following diagram commutes: 















A"+iDiff(j/) Difri_fe(i/) Difr,(^) 



Difrj_fe_i(i/) 



A* 



Difri_i(7r) 







91 







£U > 



Ker(7rf,_i)* 



The first row is a part of the Buchsbaum-Rim complex and the operator 
Vj in the second one is a differential operator with the symbol e. 

Proposition 10. Let the system £ he generalized complete intersection. Its 

compatibility on the level of {i — 1) jets, i.e. Kcr(7rf = 0, is equivalent to 
existence of an operator Vj , such that the second row is a complex. 

Proof. By the assumption the first row is exact. The diagram chase, not 
involving , yields a homomorphism 

C, : A^+^ST^i^VKeiif-Ime — > Kcr(7rf;,_i)*. (11) 

It is always an epimorphism, but it is a monomorphism iff the map A* o is 
a monomorphism (the precise form of V* is inessential here, we may consider 

its construction to be inductive by t or refer to our variety of multi-brackets). 
Vanishing of our homomorphism is equivalent to vanishing of A* o Vj . □ 

Since £.0 A* o V^* = 0, the map A* o V^* : A™+i Difr(!/, 1) -> Difri(7r, 1) from 
the proof can be identified as follows: 

j-i o A* o vr : A™+i Difr(z., 1) ^ Difri_i(7r, 1), 

which can be varied by a map of the form 

A* o Vr_i : A-+,i Diff(i., 1) ^ Difr,_i(7r, 1). 
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Since every differential operator from Diff(7r, 1) determining our system (i.e. 
vanishing on £) can be represented by operators from Diff (i/, 1) composed with 
the operator A*, the above operator represents our multi-bracket (indeed one 
of them due to non-uniqueness of Vj), or more exactly by the collection of 
reduced multi-brackets [Fi^ , . . . , = {Fi^ , . . . , F^^^^ } mod J'l_i (with an 

appropriate index i). 

The condition Ker(7rf = Vi is equivalent to vanishing of the maps 

o A* o on Aj^'*'^ Diff (z^, 1) for all t. Thus formal integrability is equivalent 
to vanishing of all multi-brackets of (m -|- l)-tuples of differential operators 
Fi, . . . ,Fr due to the system £ (as explained in theorem [A|). 

Now this line of arguments works well for the case of non-linear vector dif- 
ferential operators F — {Fi, . . . ,Fr): 

Proof of theorem |Al The above diagram over M should be considered 
over the equation £ and the Diff(l, l)-module Diff(7r, 1) of linear differential 
operators should be changed to the ^Diff'*^(l, l)-module '^Diff'^(7r, 1) of 
differential operators with coefficients in (see §2.4p . 

From the arising diagram of non-linear complexes we also obtain the multi- 
brackets as the obstructions to formal integrability and observe that their van- 
ishing due to the system ensures this integrability. 

Namely when the conditions of definition [1] are satisfied, we obtain the fol- 
lowing sequence of non- linear differential operators, which is exact at the terms 
'^Diff'^(7r, 1) and £* (in the sense specified in g231): 

A™+VDiff^(z^,l) ^'^Diff^(!^,l) -^'rDiff^(7r,l) — > £* ^ 0. (12) 

The composition £p o V* corresponds to multi-brackets of non- linear differ- 
ential operators and this sequence is a complex iff the brackets are zero due to 
the system. Thus vanishing of the multi-brackets yields formal integrability. 

Let us give more details. Compatibility on the level of i-jets, which is equiv- 
alent to injectivity of the map 7r*j_]^ : S^f„i ^ 5^f , can be expressed by saying 
that no relation of functions from Sfi{Fi, . . . , Fr) over is an operator from 
diffi_i(7r,l), see ©. 

Linearization reduces this to the following claim (cf. Proposition [9]): 

^Difff (tt, 1) n '^Diff,_i(7r, 1) = <^Difff_i(7r, 1). 

Thus due to exact sequences ([5]) and ^ the formal integrability can be expressed 
via 1-syzygy of the module £* , which is given by sequence (fT2|) . 

Indeed any non-linear relation can be evaluated at points of the equation 
Xi G £i, which gives a relation for the symbolic system g{xi). Since all such 
relations are described via the Buchsbaum-Rim complex, we describe 1-syzygy 
of £* in exactly the same way as in §4.11 

This yields multi-brackets [Fi^ , • ■ • , ^i„,+i]£ as obstructions to compatibility 
and theorem [X] is proved. □ 

A generalization. What happens to the compatibility criterion, when the 
system is not of generalized complete intersection type? 
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Then the first row of the diagram from N4.3I needs not to be exact and 
homomorphism (fTTj) does not tell much. In proposition [10] we assumed the 
system is of generalized complete intersection type, which gives Ime = Kenp*. 

If this does not hold, we get a cohomology group H^z'' at the term S^~''T®v* 
of the Buchsbaum-Rim complex . Then homomorphism ()11|) will be changed 
to 

l^ = C^® C ■■ Ker^* = Ime ® Hi,-'' — > Ker(7r|;,_i)*. 

The first component Q is represented by multi-brackets of differential operators 
as before, but the second component is of completely different nature. Again 
the map (^i is an epimorphism, but neither Q nor (^'^ needs to be a monomorphism. 

Nevertheless, since the multi-brackets together with ImC' span Ker(7r,f 
they give all compatibility conditions of the system £. 

The image of the projection TTf^_i : £i — + £i-i is the locus of the Weyl tensors 
Wi-i{S) e (we wiU not describe them here, see [Li]; for geometric 

structures they are also called structural functions [St], curvature tensors and 
sometimes torsion tensors [BCGf]). We let W{£) = ®Wt{£). Note that the 
relation W{£) = is precisely the set of all compatibility conditions. 

Epimorphism leads to the decomposition 

and so we get: 

Theorem 11. The Weyl tensor can be decomposed W{£) — Wb{£) + Wh{£)- 
The bracket part WBi£) can be expressed via multi-brackets by the same formula 
as in corollary\^ The second part Wh{£) is a homological term. □ 



4.4. An exact sequence of non-linear differential operators 

We would like now to make the complex from the proof of theorem |^ into 
an exact sequence by interchanging the functors and Diff: 

Theorem 12. Assume that the system £ is oj generalized complete intersection 
type and the multi-brackets vanish as in theorem [21 Then there is an exact 
complex of differential operators, which is a resolution of the module £* : 

^ Diff ^ (g) 5^ vr*, 1) Diff ^(A™+2zy(g)7r*, 1) 

^ ■^Diff^(A™+V, 1) ^ '^#'Diff^(i/, 1) '^DiS^iTT, 1) — >£* ^ 0. (13) 

Proof. Consider at first the linear case, where we shall use the notation 
A = {Fi, . . . , Fr) for a given vector- valued operator, defining £. 
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We have the following diagram with the exact columns: 




93 i 



<?li|pi 







Difft(A'"+iz^) ) Diff,_fe(j/) BiSiin) °' £* > 



.Difft_i(A'"+ii/) Diff,_fc_i(z/) Diff,_i(7r) £*_^ 



















The first row is an exact sequence. We will now inductively construct the 
operators Vj , □* etc, so that all the rows are exact and the diagram commutes. 

Let us choose a connection on the base manifold M and on the bundles tt 
and v. Then we can canonically embed the symbol bundle into the bundle of 
differential operators ([B [La]) via some maps g^, whence the splitting of the 
middle row of the diagram into direct sum of the first and the third rows. 

Now elements of the middle complex are 2- vectors, with the upper and lower 
components being elements of the first and the third complex respectively, so 
that the maps and A* must have the following matrix form 



e 



A* = 







where the map Vi : ® v* Diffj_i(7r, 1) is given by the system and the 

splittings and the map Xt : K'^-^^ST ® v* Diff,i_fc_i 1) is unknown. 

By inductive assumptions, the equality A* o Vj = is equivalent to the 
following: 

r,£ + A*_iXt = o. 

Thus existence of Xt is equivalent to the fact that Vte belongs to the image of 
A*_]^. By the inductive assumption this in turn is equivalent to the equality 
ai-iTiE = 0. 

The map Ti is given by the condition jiTi ~ A*q2 ~ Thus we get the 

required equality: 

{T^li-iYc^i-i^i^ = oiijiT.e = aA*q2£ - OLiqiip*e = 

because the map {'Kf^_l)* is a monomorphism and the the second row is exact 
on the level of the right three arrows. Therefore we constructed the map Vj in 
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the second row, so that the constructed part (right four arrows) is a complex in 
the commutative diagram. 

Now by the same procedure we construct successively the arrows etc 
with the condition that the second row is a complex and the diagram commutes. 

Since the first and the third complexes are exact, the diagram chase yields 
exactness of the middle complex. This proves the theorem in the linear case. 

In the non-linear case the proof is basically the same, but we need to work 
over the infinitely prolonged equation c J°°(7r), use linearizations of dif- 
ferential operators and change the rings as in ^ ^2.41 Then we get the complex 
similar to the above, evaluation of which at the corresponding jet equals the 
complex for the case of the linearized system. □ 

We will now turn the exact sequence from the theorem into compatibility 
complex for the differential operator A. For the linear differential operators it is 
known [Sl IKLV| IV] that the cohomology of the compatibility complex coincide 
with the cohomology of the stable Spencer D-complex. We can generalize this 
now to the non-linear case, whenever the system is of generalized complete 
intersection type and is formally integrable. 

Remark 10. In fact, the general non-linear Spencer D-complex is unknown, 
but in the considered situation we can construct it via the splitting as in the 
proof above. Indeed, the symbolic part is given by the diagram from the following 
section lKT] and then the complexes can be constructed inductively in the filtration 
number by the technique of the above proof. 

4.5. The compatibility complex 

The homomorphisms of the spaces of differential operators as in theorem [T^] 
are not always coming from the actual differential operators (as we treated in 
but in our situation we can arrange this. 

Indeed, in the proof from the preceding section we need to construct the 
operators Vj,^^;,^^ etc only when they appear for the first time (the corre- 
sponding index is zero or equivalently, when the defining differential operators 
are introduced). Afterwards, we can let V* = v|i\ to be the prolongation, 

— O^il-i and so on. 

Thus we get actual differential operators between vector bundles, forming 
compatibility complex for the differential operator A (or the system £). 

Let us consider at first the case of linear system £. Denote by Ss = Sols the 
sheaf of local solutions of the PDE system £. Then by the preceding discussion 
we constructed the complex 

O^Ss^ (tt) (i^) ^ (A-+V) ^ (^* ® A^+2,.) ->..., 

which is formally exact and thus is a compatibility complex of the PDE system 
£ (in some cases this complex is required to be constructed to both sides, but 
we will work only with right differential resolutions). 
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Remark 11. In the classical geometric theory of PDEs, the compatibility reso- 
lution is constructed for any formally integrable system of PDEs ( see the section 
about Kuranishi theorem in IBI), but then the construction is quite sophisticated 
(as in the second Spencer complex fSllKLVf ). In our case of generalized complete 
intersection the spaces in the complex are explicit and simple. 

In fact, the same construction will work for any formally integrable system, 
where we have managed to write explicitly a resolution of the symbolic module. 

What is more interesting is that we can write the whole compatibility com- 
plex explicitly. Consider at first two examples. 

1. Let us write exact form of the compatibility conditions for linear non- 
homogeneous scalar PDEs of complete intersection type (m = 1 and 1 < r < n) : 

Ai(u) = /i,...,A,(u) = /,. (14) 

We solve non-homogeneous linear system under assumption that the linear ho- 
mogeneous system is compatible. By theorem |A] this means certain commuta- 
tion relation between operators and their differential corollaries (we use the 
standard summation rule by the repeated indices): 

[A„A,] = C^.Afc, (15) 

^ [A,:,[A„Afe]] = 0^ [A,,Cf,]= ^ + 

cyclic: ijk cyclic: ijk cyclic: ijk 

where R'^^ = ApS^ - A^d"/! - C|^, etc. (16) 

Here C^j are certain scalar differential operators on the base manifold M" of 
ord C^j < k + Ij — Ik, Is ~ ord A^ , and are scalar differential operators of 
ordD^^]. < + Ij + Ik - min[/s +lt : s,t e {a, (3, j},s^t]-l etc. 

We may assume that is skew-symmetric by ij and D^^^, is skew-symmetric 
by ijk and by Pj. Similar conditions will be imposed on other differential 
operators with multi-indices, which occur in the higher Jacobi identities. 

Because of (|15p the compatibility conditions for non-linear system ([14]) are 

[A, - f,, A, - f,] - C^(Afc - fk) = A,/, - A,/, -I- C^k = 0. 

Thus the differential 1-syzygy operator is V = (Vij)i<i<j<„, where Vy = 
AiPj — AjPi — C^jPk with Pk being the projector to the fc-th component. In 
other words, for / = (/i, . . . , frY we have: 

V(/i, ...,fr) = {..., AJ, - A,f, - C^k, ■■■)■ 

To obtain differential 2-syzygy we write compatibility conditions for the 

system Vijif) = Xij, which are: 

cyclic: ijk 

= E (C^C|,-[A,,C'^.]-fD(^li?3,)P„ = 

cyclic: ijk 
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due to (fTC]) . Thus the next operator in the differential resolution is □ = 
{Dijk)i<i<j<k<n, where for A = (. . . , Ay , . . . ) we have 

cyclic: ijk 

Continuing in this way we arrive at the following compatibility complex: 

O^Ss — > C°°(l) ^ C°°(!^) ^ C°°(A2iy) ^ C°°(A3iy) . . . , (17) 

where z/ is a trivial bundle over M with the fiber M'' (naturally A°iy ~ 1). 

The operator : C°°(A'"'z/) ^ C°°(A'^+^z^) in the above complex has the 
following form. Let a; = {«i < • • • < ik}, ^ = {ji < ■ ■ ■ < jk+i} be ordered 
multi-indices, so that the elements of the domain for Q are 6 = {6^^}, while the 
elements for the range are 9 = {9?}- The map is 



3---Jfc + l j 



(Z?, . . . , K are some differential operators with multi-indices mentioned above 
and >Ci^k some real numbers calculated recursively, I < i < k; for instance the 
first numbers are xi^i = — i, — —1, ><2,2 = 1, xi.3 — — § and so on). 

Notice that on the symbolic level ([T7]) is the usual de Rham complex, but 
due to non-commutativity of the components of the operator A, defining the 
system £, we need to compensate the syzygy maps with the tails of (lower order) 
differential operators. Horizontal de Rham complex ( [KLV] ) is a particular case 
of (|17|). when the operators Ai commute. 

2. Consider now a vector system of (to -1-1) equations on m unknown 
functions, given by the operators Ai = (A^i, • ■ • , Aim), which is of generalized 
complete intersection type: 

Aii(ui) H h Ai„i(w„j) = /i 

A(m+l)l(ui) -I h A(m+i)„j(Um) = fm+1 ■ 

In this case the compatibility condition for the homogeneous system is given 
by the condition from Theorem 1X1 



{Ai,...,A,„+i} = 5'A, 



where {. . .} is the multi- bracket from ^3.2\ and some scalar differential op- 
erators on M of ord(_B') < J2j^i ord(Aj). 

Assuming formal integr ability for the homogeneous system, we get the com- 
patibility complex: 

O^Se^ C°°(M,R") C°°(M,]R™+i) ^ C°°(M,R) ^ 0, 
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where A = (Ay ) is the (m + 1) x m matrix determining the system £ and 

m+l 



V(/i 



= E (-l)'-'Ndet[A,,]J|^'^™ -SM A 



fc=i 



For m — 2 this formula looks as follows (the determinants are non-commutative): 



A21 A22 
A31 A32 



+ 



A31 A32 

All A12 



All A12 
A21 A22 



3. Now we combine the first case of arbitrary dimension of base with the 
second case of arbitrary dimension of the fiber of tt as in §4.21 Then we can write 
the compatibility complex explicitly (in the generalized complete intersection 
case; in general it is more complicated to write the W^i/(£)-component). 

Let 77 and i/ be trivial vector bundles (locally this is not a restriction) . More- 
over we will assume v is graded in the various orders case for the operator 
A : C°°(7r) -> C°°{iy). 

Theorem 13. The compatibility complex for a collection of linear operators 
A — {Ai}, defining a system £ of generalized complete intersection type, has 
the following form: 



O^Se^ (tt) ^ Cr^,{u) ^ 



Q~ (A^+i^.) ^ 



<--ioc(7r ^lA ^ V) > C;o^(A TT 



A 



m+3. 



The differential 1-syzygy operator V : C°°(M] 
by the formula 



C°°(M; A™+iR'~) is given 



m+l 



V.(/) = E(-l)'"'Ndet[A,,,]^|P"/, 

fe=i j=i 



Em 



where f — (/i, . . . , fr) and V ~ {Vr}; t being a multi-index 1 < ii < ■ ■ ■ < 
im+i < f- The coefficients B!^ are obtained from the multi-bracket relation 
(compatibility of the linear homogeneous system) 



{A,,,...,A,„^J 



The higher syzygy operators are given by the method above. 



□ 



Non-linear case is treated by the same formulas via the linearization op- 
erator. We need to restrict to infinitely prolonged equation c J°°(7r). 
Denote the projection from this equation to the base M by TTg. Then changing 
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the operator A to its linearization £a (the operator itself on the equation van- 
ishes) and coupling this to the arguments above we will arrive at the non-linear 
compatibility complex 

(note the change of the solution sheaf Ss — Solg in the linear case to the 
symmetries sheaf Symg — Ker(£A) in the non- linear case: For a linear equation 
£ each shift by a solution is a symmetry). 

This complex is formally exact and its cohomology can be identified, as 
we have noticed in the previous section, with what we call non-linear Spencer 
cohomology. 

This makes an effective representation for this important invariant H\~i{£) 
of the system £ of differential equations (recall that in general non-linear case 
this invariant is not defined). 

5. Integrability theory 

In this section we consider certain topics closely related to the main subject 
of this paper. In fact, explicit compatibility/solvability criteria can be applied 
to integrability of determined PDEs. Indeed, one can find (simple Frobenius 
or sometimes more sophisticated) compatibility schemes in most well-known 
integrability approaches, in particular in the following theories: 

— Symmetry calculus; Backlund transformations; 

— Lax pairs; Zero-curvature representations; 

— Darboux integrability; Sato theory etc. 

Thus it is important to understand formal integrability and related Poisson 
geometry, method of differential constraints via multi-brackets or other criteria. 
Some of these topics will be considered in this section. 

5.1. Spencer cohomology and curvature tensors 

Let £ C J°°(7r) be a system of differential equations defined by a set of 
differential operators A : C°°(7r) C°°{v), which we allow to be of different 
orders. 

Proof of theorem |Bj Below we denote hy g ~ ®gi the symbolic system 
associated to £. 
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Consider the following commutative diagram: 



0^ g -> ST* ®7r ST* (g) v A^+^ 

0^ g(E)T* ST* (E) TT (E) T* ST*(E)v(E)T* A™+i ® T 



0^ g(g) A^T* ST* TT A^T* ST* (g> v (g> A^T* A^+i ® A^T* 



0^ g(g) A"T* -> ST* (g)Tr(g) A"T* ^ ST* (g) v (g) A"T* A™+i (g) A"T* ^ ... 

\- \' \' \- 



where we write A™+^ = ST* g) A^+^z^ for brevity. The rows are generalized 
Spencer complexes (i.e. M-dual to Buchsbaum-Rim) and are exact. 

The columns are the usual Spencer (5-complexes and all except the first one 
are acyclic, i.e. have cohomology in the bi-degree (0, 0) only. Thus by the 
diagram chase the cohomology of the first column at the j-th term H*'^{g) 
equals the zero cohomology group of the (j + l)-th column. □ 

Denote /i*'-' ~ dimH'^'^ (g) the Betti numbers for Spencer cohomology. 

Corollary 5. For a generalized complete intersection g we have: h*'^ — m, 
= r, h*'^ = r^lf)^;^.^) for2 < J < r + l-m and h*^^ = for 
r + l — m<j<n. □ 

Note that the Euler characteristic vanishes as it should: x = J2i~^y f^*'^ — ^■ 

It is possible to specify the bi-grades, where the J-cohomology does not 
vanish. If g is a system of pure order fc, then the non-zero Betti numbers are: 
h°'°, /I'^'-i'i and for 2<j<r. 

For a system of different orders 1{1), . . . ,l{r) the above complex allows to 
specify the bi-grades too. The precise combinatorics is straightforward from the 
weighted Buchsbaum-Rim complex. However, since the formulas are involved, 
we indicate only what happens to first and second (5-cohomology (this latter 
space contains curvature, i.e. Weyl tensors, see Corollary (2). 

For them the only non-zero Betti numbers are: /i'^*)"-'^'^, 1 < i < r and 

/j/(*l) + -+/(Wi)-2,2^ 1 < il < • • • < im+1 < r. 

Remark 12. Formulas for h*'^ from corollary\B[ suggest some polynomial rela- 
tions between multi-brackets. In fact, form ~ 1 we can see the basis of H*-^ for 
J > 2 as the set of power j subsets in {1, ... ,r}: for each such subset we asso- 
ciate all possible iterated Mayer-Jacobi brackets, which is an analog of j -form. 
The Jacobi identity and its higher analogs yield relations between these iterated 
brackets {...{{F,„ F,,}, F,J, . . . , Fi^. 
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When m > 1 the number of subsets of power m + j — 1 in {1, ... ,r} is 
im+j-i)' -'Another factor o/dimiJ*'-' forj > 2 is related to the fact that iteration 
of brackets is now arranged in a multiplicative manner, so that the relations are 
generalized Plucker identities as in ^3.3[ 

5.2. Integrability of characteristics and multi-brackets 

Consider the space of scalar linear differential operators Diff(l,l) on the 
manifold M. It bears the structure of infinite-dimensional Lie algebra when 
equipped with the Jacobi bracket {, } (this is also true for the space Diff (tt, tt)). 
This bracket induces the classical Poisson bracket on T*M as follows. 

Let a : Diff(l, 1) C°° {T* AI) be the symbol map. It associates to an order 
k differential operator F a polynomial a{F) — smblfe(-F) of degree k in momenta 
p in canonical coordinates (x,p) on T*M. Then we have (with an ambiguity in 
notations): 

MF),a(G)} = a({F,G}), 

where the brackets on the left are Poisson and on the right are Jacobi brackets. 
In order to make a distinction we will write {, }^ for the symbolic (2-)bracket. 

Consider a system of linear scalar PDEs £ ~ {Fi ~ 0, . . . , F,. — 0}, Fi £ 
Diff(l, 1). Let T = T{F) be the differential ideal generated by Fi, 1 < i < r. 

Proposition 14. Let the scalar system £ be formally integrable. Then the cor- 
responding ideal X is closed with respect to the higher Jacobi bracket. 

Proof. The system £ is defined by the ideal I whenever formally integrable. 
Jacobi bracket of two operators is their differential corollary. Since the ideal T 
is differentially closed the claim follows. □ 

Corollary 6. The corresponding .symbolic (characteristic) ideal I — is 
closed with respect to Poisson bracket. □ 

A stronger statement was proved in I GQS| : Namely the radical \/T is Poisson- 
closed too (in fact, the claim was justified only for the components of Char'^(5) 
of maximal dimension). This is the celebrated integrability of characteristics: 
The affine characteristic variety Char^g(5) C T*M is integrable in the Frobe- 
nius sense, i.e. if F and G vanish on it, their Poisson bracket {F, G} vanishes 
on it as well. 

This was applied in 'GQS| to the case when T (resp. /) is the annihilator 
of the module £* defining a formally integrable PDE system (resp. the sym- 
bolic module g*). However, for the system involving several unknown functions 
the characteristic variety does not bear the complete information about the dy- 
namics in generalization of Hamilton- Jacobi theory. Indeed, as we described in 
§1.31 the characteristic variety Char^(f ) of the system £ defined by a differen- 
tial operator A : C°°(7r) C'^{v) is the support of the characteristic sheaf 
K = Kercr(A). 
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Dualizing this we can equally work with Cokera(A*). We describe at first 
the picture with differential operators. The Diff (1, l)-module £* is given by the 
exact sequence 

Diff (zy, 1) ^ Diff (tt, 1)^E* ^ 0. 

Denote J = lni(i^'^) the submodule in Diff(7r, 1). 

The following statement is obtained similarly to proposition 1141 

Proposition 15. Let the vector system £ he formally integrable. Then the cor- 
responding module J is closed with respect to the multi-bracket {•,...,•}. □ 

This is the necessary condition for formal integrability from theorem]^ It 
is sufficient only in special cases, see below. We now define the symbolic 
multi-bracket generalizing the above Poisson 2-bracket. 

Definition 4. The symbolic multi-bracket of sections fi = (^{Fi) is given by the 
formula: 

. . . , (7(F„+i)}, = (7({Fi, . . . , 
(we use the index a just to keep distinction between multi-brackets) . 

It is easy to see that the right-hand side does not depend on the lower 
order terms of Fi, so that the left-hand side is the well-defined expression 
{/i, . . . , /m+i}(T- However the above formula gives the symbolic multi-bracket 
only for vector- valued polynomials / £ ST^tt* . The standard trick extends this 
to formal series and analytic functions, but we can define the symbolic multi- 
bracket {• ■ • }(T to all smooth vector- valued functions on the cotangent bundle 
fi e C°°(r*M; TT*), 1 < i < TO -I- 1. This follows from the following proposition. 

Let Ci be a basis in the bundle tt* (which we assume trivial or make a 
localization). Then a vector-valued function / e C°°{T*M ; tt*) can be identified 
with the collection of functions € C°° {T* M) via the decomposition / = f^ej. 
We use the components below. 

Proposition 16. The symbolic multi-bracket of the vector-valued functions fi 
can be expressed via the product and the standard Poisson bracket of their com- 
ponents ff (1 < i < m + 1 . 1 < i < m) . 

Proof. The multi-brackets of differential operators Fi = {F^, . . . , FJ^) has 
the following form in components: 

{Fi, . . . = ^ (-1)^F^^(i)F^(2) • ■ ■ K{7n}K(m+l)- 

Taking the symbol and making elementary transformations we get (both indices 
j and s vary between 1 and to) the multi-brackets of fi — {fl, . . . , /™): 

{/l; • ■ • J /m-l-l}^ = (-1)™ ' ^ [2 n-^^(«) ' ''-•^^(0''^T-(m+l)}cr 

rGSm+i s^l 

+ X] n -^Tis) ■ {/!(()' Ku)}'^ ' /r(m+l) (18) 
3>l s=i],l 
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This represents components of the symbolic multi-bracket via the 2-bracket. □ 

Remark that our symbolic multi-bracket differs from other multi-versions of 
the Poisson bracket, like Fillipov-Nambu or generalized Poisson bracket ( |APB|, 
IMVl [N] ). For instance, the multi-bracket with m > 1 is not a derivation in its 
arguments. However it is not pretty far from this: 

Corollary 7. The symbolic multi-bracket {/i, . . . , fm+i}a- is a differential op- 
erator of the first order in each of its arguments fi G C°° {T* M , tt*). □ 

Indeed, we can write for g g C°°(r*M) 

{/17 • ■ • J gfii •■•II /m.+ l}(T — 3{/li • • • I /ii • • • I I fm+l}a 

= E,, fl}a ■ •••,/.,.•., I fm+l), (19) 

where the exact form of the functions c^^ can be obtained from (|18p . 

Let J — (7{J) be the symbol of the submodule J . Denoting as in i i2.2ll2.3l 
= (TA the symbol of the differential operator A, we can define J via the exact 
sequence 

ST ®v* ^ ST ®iT* ^ g* -^Q 

as J = Iu\{lp) C ST ® TT*, so that the symbohc module is g* = ST ®ti*/J. 
Then proposition 1151 vields: 

Corollary 8. The (characteristic) submodule J = cr(j7) is closed with respect 
to the symbolic multi-bracket. □ 

Beside identity the symbolic multi-bracket satisfies the same properties 
as the multi-bracket of differential operators, for instance we have the general- 
ized Pliicker identity as a corollary of ^3.31 (the upper index means component): 

'^^{~^)^{{h, • • • I /fci • • • I fm+lVay ■ fk = ^'^{—^)''{fl, • • • I A, • • • I /m+2}(T ' /fc, 

where {• • • }J is the symbolic multi-bracket associated to the opposite multi- 
bracket of differential operators {• • • }^ of i i3.2l bv the same rule as in definitionlH 

Consider again the case (of importance to PDEs), when fi <E ST (g) tt* are 
polynomial vector-valued functions on T* M . 

Let /o(J) C i? = ST be the ideal generated by the functions // and Co its 
Poisson center. The following statement follows directly from (fT9)) : 

Corollary 9. The multi-bracket is a homomorphism in its arguments over the 
ideal Co dST: {■■■}„ ^ Homco (A™+i J, J) . □ 

The symbolic multi-bracket can be used to formulate a Hamiltonian 

formalism as generalized Poisson brackets are used. Namely, let /i, ... , fm € J 
be Hamiltonians and h E ST (g) a* a polynomial vector-function. A multi- 
Hamiltonian operator 

^/iA...A/„. : g* g*, h mod J 1-^ {h, /i, . . . , f„i}cr mod J 
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determines transport on Char (g) (symbolic module g* is supported on char- 
acteristics). Caustics of solutions (wave fronts) develop according to it. More 
details will be provided elsewhere. 

5.3. Applications to smooth integrability of PDEs 

The classical Lagrange-Charpit method [Gou| IGuj is designed for first order 
scalar PDEs and is as follows. Let F = Fi{x^ , . . . ,x''\u,pi, . . . ,pn) = be 
a differential equation. To solve it one searches for functions F2 , . . . , F„ on 
J^(R") such that [Fi,Fj] = mod (Fi, . . . ,Fn), where the [, ] is the classical 
Mayer bracket. Then the system Fi = 0, . . . , F„ = gives a finite dimensional 
family of solutions of the PDE F = (this family is equivalent to an ODE). 

If in addition F2,...,F„ are symmetries of the system, [F,Fi\ = X ■ F, 
then one can obtain a complete integral of the PDE F = as the system 
Fi = 0, F2 = C2, ■ ■ ■ , Fn = Cn, which in addition can be found in quadratures by 
the symmetry method of S.Lie ( jLiej ). 

Example 1. Consider the PDE upi ■ ■ -pn = ■ ■ ■ x" on R". It possesses a 

collection 0/ auxiliary integrals; F,; — — — — . This gives a complete integral of 
the differential equation. 

Now basing on our theorem \^ we can formulate a generalized Lagrange- 
Charpit method. For manifolds of dimension 2 this was done in IKL2I . To 
obtain the general version we start with the following idea. 

Definition 5. Let £ be a formally integrable system of PDEs. Call a system £ 
an auxiliary integral (or a set of integrals ) for the system £ if the joint system 
£ n £ is also formally integrable (= compatible). 

We proved in IKL2I that classical objects, such as point symmetries, contact 
symmetries and intermediate integrals f [Gou[ [Gul ILE| ) as well as higher sym- 
metries { |KLV| ) are partial cases of this notion. Moreover, some of the newly 
introduced generalized symmetries are also auxiliary integrals. 

Another traditional method for finding exact solutions of PDEs is a method 
of differential constraints. Then one considers an overdetermination f on a sys- 
tem £, such that the system £ C]£ is solvable. The solvability is very non-trivial 
to check in practice. Thus our notion of auxiliary integral is more constructive, 
since we can use an effective criterion from theorem |X] to check compatibility 
(see |KLi[ IKL2I for examples) . 

Generalized Lagrange-Charpit method is the following special form of an 
auxiliary integral. Consider a determined system of m PDEs £ ~ {Fi = 
0, . . . , Fya = 0} on TO unknowns functions ui, . . . , (we can even start with 
an underdetermined system). We search to add to it n — 1 differential equa- 
tions £ = {F„j+i = 0, . . . , F„i-|_„_i — 0}, so that the resulting overdetermined 
systems £ {^£ = {Fi = 0, . . . , F„i+„_i = 0} is: 

a) generalized complete intersection; 
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b) compatible. 

Of course, one can add less functions, but advantage of (n— 1) is that the system 
becomes of finite type (being compatible it constitutes an integrable distribution 
by the Frobenius theorem St, KLa l) and thus reduces to a system of ODEs. 
By theorem |A] compatibility of £ n £ is given by the conditions 

for some differential operators A of orders ord(A^^ J < + - ■ ■+im+i—ij~l. 

On the symbolic level the first step is to include / = o'(F) into a submodule 
J — (/i, /2, . . . , /m+n-i) of ST ® TT*, which is closed under the symboHc multi- 
bracket, as integrability of characteristics from ^ 35. 21 claims: 

{/in ■ ■ • I /im+i}(T = ■ fj- 

Then we shall adjust sub-principal symbols of Fi. 

Example 2. Consider the Cauchy-Riemann system on the plane: 

£ : = Vy, Uy = -Vx- (20) 
Then the differential equation 



£ : det 



Ux Vx 



G{u,v) 



^y Uy^ 

is an auxiliary integral suitable for generalized Lagrange- Charpit method iff 

AG = ||VG||VG, (21) 

where A is the standard Laplacian and V is the standard gradient on the Eu- 
clidean plane M.'^{u,v). Thus our transformation can be seen as a kind of 
Backliind transformation, which for any solution G{u, v) of equation \21\} as- 
sociates a 3-dimensional family of solutions of system I120\) (for dimensional 
calculus see the next section). 



5.4. Formal dimension of the solutions space 

Consider a symboHc system g = {gi C S'T* (g) N} and let V* C T* be a 
subspace. Then we can define another symbolic system g = {giH S'-V* N} C 
5"^* (g) N. It is called the V*-reduction of g ([KLj). Denote W* = T*/V*. 

We take g to be the symbolic system of a generalized complete intersection 
£ of formal codimension r, the same number as in definition [TJ Note that 
this definition can be reformulated for symbolic systems as well. Then the 
characteristic variety Char'^(g) = Char'''(£) has codimension r — m+l in p^T*. 
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Theorem 17. Let g be a symbolic system of generalized complete intersection 
type and the subspace V* C T* of dimension (r — m + 1) be transversal to the 
affine characteristic variety of g: Cha,r{g) n P^V* = 0. Then the reduced sym- 
bolic system g is also a generalized complete intersection, Spencer S-cohomology 
of the system g and of its V* -reduction g are isomorphic and g ~ g (g) SW* . 

Proof. This statement was proved for more general Cohen-Macaulay sym- 
bolic systems in jKLgj- Since we proved in §2.31 that a generalized complete in- 
tersection is Cohen-Macaulay, the claim follows (the part that the y*-reduction 
form a generalized complete intersection is straightforward). The last claim of 
the theorem gi ~ J^j 9i-j ® S^W* is not a part of theorem A from IKL2I , but 
is contained in the proof, see remark 8 loc.cit. □ 

Notice that due to Noether normalization lemma [E] a generic subspace 
V* C T* is transversal to the characteristic variety over C. The functions on 
W — Ann(y*) are those on which a general solution oi £ depends, as is refereed 
to in the discussion of functional dimension before Theorem O 

Proof of theorem[Cj Let us consider at first the case, when r = n+m— 1. 
Then the system £ is of finite type, namely iVi^i-i : £i ^ £i-i for i > X]I=i ^i- 
The equation f 00 is thus a finite-dimensional manifold equipped with the Cartan 
distribution Ce ( |KLV|, [KLTI ) ■ Since the system is compatible, the local solutions 
of £ are integral manifolds of the distribution Cg of dimension n. 

Thus the dimension of the solutions space is 



00 



dim Sg = dim £00 — n = dim 



(the sum is indeed finite) . We can calculate dimensions of the symbol spaces gi 
explicitly, since the system is of generalized complete intersection type. 

For brevity sake we omit the immense combinatorics and provide details of 
this step only for the case n = dimM = 2. The general case is similar (we refer 
to IKL5I , where arbitrary systems of pure order k are scrutinized) . 

If n = 2, we get r — m -\- 1. Let ki < ■ ■ ■ < kr he the orders of the system g. 
We have: 



dim gi 



m{i -\- 1) if i < fci 

(m - j){i + 1) + Y^i^i ks if kj <i < kj+i 

J2l=i ks-l-i \ikr <i < X]I=i ks 

else. 
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Thus 

fel-l fe2-l 



^dim.9, ^ ^m{i + l)+Y^ (^{m - + 1) + fci) 

i—O i—ki 
ka-l k,.~l 
+ ^ ((m - 2){i + 1) + (fci + fca)) + • • • + (^1 + ■ • • + >^r-l I 



i—k2 i—kr—i 

k^+■■^^-l , , , A (fci + ■ ■ ■ + kr-i - l)(fci + ■ ■ ■ + fc^^i) 

+ 2^ (fci H hkr-l-ij = 

i—kr 

+ J2 ^'^V + ^^^^ " ^'-^^^^^ + ■ ■ ■ + ^'-^^ ^ ^ 

2—1 i—1 i<j 

and the result follows. 

Consider now the case m<r<n + m— 1. Then the characteristic variety 
Char'''(f ) is non-empty and the system £ is of infinite type. By theorem [T7l the 
symbolic system g has a free factor 5*11^* and the dimension 

p = dim W ^ n — dim V = n + m — r — 1 

is clearly the formal functional dimension of Solg . 

Thus the quantity X]i=o 9i when t oo grows as d ■ dim S*W* , where d is 
the formal functional rank for the system g, which is of finite type (no complex 
characteristics) . The reduction g is of generalized complete intersection type and 
it has the same orders ki < ■ ■ ■ < kr because the first Spencer 5-cohomology 
coincide (theorem [TTl) . Thus we can use the calculations above to conclude 

d — dim cji = ki^ - ■ ■ ki^ — d, I = dim V = r — m + 1. 

ii<---<ii 

Since dim W* = n + m — r—l, the asymptotic of the Hilbert polynomial for the 
symbolic module g* is 

d-dim5*M^* ' 



(n + m — r — 2)! 

This proves the theorem. □ 

Corollary 10. // a generalized complete intersection £ oj formal codimension 
r is formally integrable and has equations of the same order k, then its for- 
mal functional dimension and functional rank are p = n + m — r ~ 1 and 
d = ("'''™~^)^" (i-S- o, general formal solution depends on d functions of p 
arguments). □ 
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5.5. Generalizations of compatibility via brackets 

Vanishing of multi-brackets is necessary, but not a sufficient condition for 
compatibility. By theorem [A] it is sufficient for generalized complete intersec- 
tions, however this does not generalize to more general class of Cohen-Macaulay 
systems. 

To see it consider a system of finite type. It is a Cohen-Macaulay system. 
In fact, the finite type condition means Char ''(g) — 0, so that dimg* — 0. But 
this condition also implies depth = (alternatively < depth g* < dimg*). 

Now for a completely determined system of order k (Frobenius type), i.e. 
gk — 0, the compatibility conditions do not coincide with these of our theorem. 
Actually, let us write the equation in the orthonomic form: 

S = {^ = F^{x,^)\l<^,J<m, \a\^k,0<\r\<k}. 

The compatibihty condition is <i>^ f,^^ = 'DaF^_^^^ - VbFi^^^ = mod Jk{S)^ 
\t\ — k — 1, while the multi-brackets are: 

- , . . . - Fl'^^l\} = Ht + (smaller order terms), 

where the higher order term Hx (with order equal to I'^'il ~ 1) non-zero 

iff for some indices a, (3 we have (in the display below the braces mean a set): 

«Q = i/3 = fc e {1, • ■ • , "} and {ii, . . . ,ia, . . . ,if3, . . . , im+i} = {1, . . . , fc, . . . , n}, 

in which case = Sl{±Ilj^a,p'Daj)['Da^F^^ — V^i^F^J. This provides more 
conditions and they are of higher order. 

Remark 13. It is possible however to give explicit compatibility conditions for 
some systems different from generalized complete intersection. Let £ be a Cohen- 
Macaulay system and g its symbolic module. Choosing a subspace V* C T* not 
meeting the characteristic variety Q}l^x^{g) and of complimentary dimension, 
we get the reduction g, which is also Cohen-Macaulay by theorem A of JKL^. 

If we have a compatibility criterion for systems £ of type g, we can transform 
it to get a criterion for the system £ of type g (cf. theorem \17\ l. For instance, 
if the system g is completely determined (Frobenius type), we can use the above 
formulas to get similar compatibility conditions for the system £ (which is not 
of Frobenius type!). 

Thus usually for systems different from generalized complete intersections, 
multi-brackets do not provide a basis of compatibility conditions (though multi- 
brackets are part of them) . Indeed, in these other cases the obstructions to for- 
mal integrability (Weyl tensors) belong to different Spencer cohomology groups. 

However in some cases the system not of generalized complete intersection 
type can have compatibility conditions in a form of multi-brackets. Usually this 
happens when the Spencer cohomology is of the type described in theorem [b1 
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For instance, for a (skew-)product of a generalized complete intersection and 
some involutive system. We give two examples. 

1. Let J be an almost complex structure on a manifold M. It defines the 
Cauchy-Riemann operator dj : ^^''^{M) — > il.P'^{M), so that the system £ = 
Ker{dj) is not a generalized complete intersection. However we can represent 
the CR-operator locally as the product of the operator Bj : C°°{M) — > ^W'^IM) 
and the identity on 17^'° (M). The compatibility condition (which is equivalent 
to integrability of the structure J) has now the form of vanishing multi-brackets. 

2. Let V be a connection on the bundle tt : E ^ M and ^^^•'^{E) the 
bundle of p- vertical, q-horizontal forms. The horizontal de Rham operator : 
nP'°{E) nP'i(£') is locally a product of nP^°{E) and the horizontal de Rham 
operator dy '■ C°°{tt) ilP'^{E). The latter has compatibility condition (which 
is equivalent to flatness of the connection V) again in the form of vanishing 
multi-brackets as in theorem |Al 

6. Applications 

In this section we apply the compatibility criterion to solve certain problems 
arising in differential geometry. Many of them address solvability of overdeter- 
mined systems of PDEs. 

To decide if the system is solvable we add compatibility conditions to the 
system (here we use theorem]^, investigate the new system, add its compatibil- 
ity conditions etc. In other words we apply prolongation-projection scheme and 
either close up the system or get a contradiction (empty equation £ C J'^{tt)). 

For example, linearization problem of 3-webs on the plane is equivalent to 
solvability of a system of 2 scalar second order equations of two variables. The 
system is a complete intersection and the compatibility condition is given by 
vanishing of Jacobi-Mayer bracket, as was sketched in |KLi| . This condition is 
equivalent to vanishing of the Chern curvature and yields parallelizable webs. 

To linearize the web one adds the bracket to the system and further in- 
vestigates compatibilities. This was done in |GLi| and thus the long standing 
Blaschke problem was solved. 

Another example is the problem of finding the number of Abelian relations, 
which is equivalent to solvability of a system of (to + 1) differential equations 
on TO unknown functions on the plane. The system is a generalized complete 
intersection and the compatibility conditions is given by vanishing of multi- 
brackets. This method was applied in IGL2I and the rank problem, addressed 
by Lie, Poincare and Bol, was solved. 

In this section we solve with our technique some other pending problems of 
classical differential geometry. 

6.1. Killing vector fields on the plane 

If a Riemannian metric 5 on a surface possesses a Killing vector field, it 
has the following local form: ds^ = gii{x)dx^ + 2gi2{x)dxdy -\- g22{x)dy^ (near 
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the point, where the field does not vanish) and vise versa, so that this is a 
surface of revolution. 

Now we address the following question: How to recognize if a metric can be 
brought to such a form? This classical question was studied by Darboux in [D] . 
We however did not find a clearly formulated answer in the literature. Here we 
give a criterion in differential invariants using our compatibility technique. 

The problem is equivalent to solvability of the equation L^ds^ = 0, where 
^ = udx+vdy is the required vector field and ds^ = gij{x, y)dx'^dx^ is the metric, 
a;-'^ = x,x'^ — y. The tensor equation is equivalent to the following 3 first order 
linear PDEs on the functions u{x,y),v{x,y): 

'2uxgii + 2vxgi2+u{gii)x+v{gii)y = 0,2uygi2 + 2vyg22 + u{g22)x+v{g22)y = 0, 



Denote them by Ei,E^,E2 respectively. We get the linear system E = Ei C 
J^(2;2) of codimension 3 (we shall write J^{n,m) instead of J'^(E",R™)), so 
that dimgi = 1, dim 32 = 0, whence the isomorphism 1:2,1 '■ E2 — > Ei. 

The compatibility condition is equivalent to the Frobenius condition on the 
corresponding distribution L(7r^i[) on Ei and is given by the condition E'^ = 
[Ei,E2,Es]£ = £ El- This differential operator E'^ has order 2, but due to 
the above isomorphism can be considered as a function on J^(2;2). However 
if we consider it modulo Ei — 0, E2 = 0, E3 — 0, it becomes a function on 
J'^(2; 2) (this is not automatical and is a peculiarity of the system) and has the 
form: 



where \g\ = 511522 — ffi2 is the determinant of the metric and K is the Gaussian 
curvature. Thus compatibility condition is equivalent to the claim that {M'^,g) 
is a spacial form: K = const. Note that this is the case, when the solutions 
space has dimension 3. 

Let us study solvability, then we need to add the equation = to the 
system. This means u = KyW, v — —K^w and we obtain the following system 
on one function w{x,y): 



where a = guKy - gi2Kx, (3 = guKy - g22Kx, 71 = {gii)xKy - {gxi)yKx + 
2giiKxy - 2gi2Kxx, 72 = {g22)xKy - (522)2/^1 + 2gi2Kyy - 2g22Kxy, 73 = 
{gi2)xKy — {gi2)yKx + giiKyy — g22Kxx- Note that a, /3 do not vanish simulta- 
neously unless Kx = Ky = 0. 

Denoting by Si the determinant of the above matrix we obtain two necessary 
and sufficient conditions for non-trivial solvability (uj = is always a solution): 
5*1 = and S2 = 0, where: 



Uygii + Uxgi2 + Vygi2 + Vxg22 + "(312)0; + v{gi2)y = 0. 



Ei^ E'^{modEi,E2,Ez) = ^g\''{KxU^ Kyv), 




= 



If a/3 7^ 0, then ^2 = (0(71) 



V 



ay-ii)(i^-(P{^2)x-M2)a\ 
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- If a = 0, /? ^ 0, ^2 = Pil2h - M2 - P{l3)y + f3yj3 (^1 = ^ 7i = 0). 

- If a 7^ 0, /? = 0, ^2 a(7i)a ~ "a7i - "(ts)^ + ax73 (5'i = ^ 72 = 0). 

Thus the criterion for existence of Killing vector field becomes the following 
two non-linear differential relations S'l = and 6*2 = 0, having orders 4 and 5 
in the coefficients of the metric g respectively. 

For a tensor T denote d^^T = dy^dyT) the covariant derivative of the tensor 
dyT differs fr om dy, which is equal to multiplication by the curvature 

tensor). In particular, we obtain the forms d^'if S C°°{& T*M). Note that the 
form is symmetric, but the higher covariant derivatives lack this property. 

Let also grad K be the g-gradient of the curvature and sgrad K = J grad K 
be its rotation by tt/2 (fix orientation). The preceding calculations imply: 

Theorem 18. The space of local Killing vector fields can have dimension 3, 1 
or zero. A Riemannian metric g possesses a local Killing vector field iff 

d^^K {grad K,sgTadK) = and (sgrad sgrad if, sgrad iiT) = 0. 

There are 3 independent Killing fields iff K is constant. 

Remark 14 The main claim of the theorem is the formula ( sufficiency of which 
is obvious). Other statements were known to Darhoux JD^ . In fact, even formu- 
las can he attributed to him, though no precise statement was made in fDj; see 
\K2^ for details. 

Note also that global implications are straightforward, but the dimension of 
the space of Killing vector fields can differ. For instance, for the standard flat 
torus it is 2. 

Proof. Let us note that ii K ^ const, then the system 

£' = {El =E2 = E3 = E4 = 0} 

has symbol dimensions: dimf^g = l,dim55. — and so in the compatible case 
the solution space is one-dimensional. Thus we need only to prove the existence 
part of the theorem. 

For this we express the above S'l and iS'2 via differential invariants. Direct 
calculation shows: 

S'l = -A\g\^d'^^K {grad K,sgTa.dK) and 

S2 = Adf KisgvadK, sgrad K, sgrad K) + ^^f Di\g\^Si), 

I grad K \ * 

where in isothermal coordinates, when ds^ — e^{dx'^ + dy^), we have: 

-2K^K^\g\^ , , 

A- Ig^^d^P ^"'^ U^-K^Kp^-KtKyDy+ 

+ {K^yK^,, + KtKyy + 2K,Ky{Kl + Kl)K„j - 2KIKI{KK,, -t- XyKy)). 
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Thus l^l"^^! is a differential invariant, while A~^S2 is a differential invariant 
relative the condition Si = 0. □ 

In the next section we'll need to enumerate differential invariants of a Rie- 
mannian metric on a surface. It is a known fact (see [T]), that the space of 
scalar differential invariants of order fc of a Riemannian metric on a surface is 
generated by (fc — 1) differential invariants for all fc > except fc = 3, where 
there is only one invariant. 

The first invariants are: I2 = K and I3 — \WK\'^ (the index refers to the 
order of differential invariant). 

To fix a basis in invariants of order i = 2 + fc we consider the form d^^K and 
substitute gradX as first (i — j) arguments and sgradii' as the next j arguments 
(0 < J < i). We denote the resulting function lij and enumerate the index j by 
letters (so we write I^i, instead of /41, I^d instead of 153 etc). 

In these invariants the criterion of Theorem 1181 writes : I^b — 0, Isd — 0. 

6.2. Higher order integrals of plane metrics 

Killing vector field on a surface M can be represented as a linear (in mo- 
menta) integral of the geodesic flow on T*M. It is important to know when 
the flow admits a polynomial integrals. Locally geodesic flows are integrable, 
but the corresponding integrals are usually analytic only on T* M \ M . So in 
general polynomial integrability requires certain conditions even locally (here 
and throughout the standard regularity assumption should be imposed). 

Let {x,y) be local coordinates on and Px,Py be the corresponding mo- 
menta on T*M . Since every homogeneous term of the integral is obviously an 
integral, we consider a function Fd — J2i+j=d ^iji^^ y)PxPy of degree d on T*M 
{i in p\ is a power, not index). 

The Hamiltonian of the geodesic flow is H = g^^p^+2g^^pxPy+g^^Py (matrix 
g^^ is inverse to the matrix gij of the metric). Let {H, Fd} be the Poisson bracket 
of H and Fd- It is a polynomial in momenta of degree d -I- 1. 

Thus involutivity condition {H, Fd} = is equivalent to {d + 2) equations 
£^1=0,..., Ed+2 = on (d-|-l) unknown function ado{x, y), . . . , aQd{x, y). These 
equations form the first order system E of generalized complete intersection type 
and so the compatibility condition can be expressed via the multi-bracket 

Ed+3 — [El, . . . , Ed+2]£ — 0. 

If this condition is not satisfied we add Ed+3 to the system and continue with 
investigation of solvability. 

In this section we consider the case d= 2. This is the classical case, studied 
since Darboux. It is known ([D | IKoj . see also [Bi]) that existence of an addi- 
tional integral, quadratic in momenta, is locally equivalent to the possibility of 
transforming the metric to the Liouville form 

ds^ = {f{x) + h{y)){dx^+dy^). 
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However, no effective criterion for recognizing Liouville metric was obtained 
despite many attempts. The only visible success was a note |Su| . The solution 
to the problem was sketched there, but the answer was not written in invariant 
terms (notwithstanding the title), and the number of differential invariants char- 
acterizing Liouville surfaces was not given (in fact, it is difficult to pursue what 
the proposed set of compatibility conditions actually is and why it is complete, 
so that we choose another approach below) . 

We describe a criterion basing on our compatibility criterion. Let us write the 
metric in isothermal coordinates: ds^ = e'^^^'^'(da;^ -I- dy^) (the approach works 
with the general form as well, but the expressions become too complicated; 
However since the answer will be given in differential invariants, the method 
plays no role). 

The function / — u{x, y)dx'^ + 2v{x, y)dxdy -\- w{x, y)dy^ is a quadratic inte- 
gral of the geodesic flow iff the following system £ (coefficients of {H, 1} = 0) 
is satisfied: 



Denoting the equations by Ei, E2, E^, E4 we obtain their compatibility condition 



From the general theory it might be expected that E5 has order 2 (3 in non- 
reduced form), but in fact it is of the first order and has the following form 
(after cancelation by 2^ydet{g)): 

E5 = 5KxVx - 5KyVy 

- (Kxy + 2\yKx + 2\xKy){u - w) + {Kxx - Kyy + AX^ K ^ - ■iXyKy)v. 

Thus the system £ is integrable iS K ^ const. In this case dimension of 
the solutions space is ^dimgk = 6 {dimgi — max{3 — i,0}) and the space of 
quadratic integrals is the symmetric square of the 3-dimensional space of linear 
integrals (a basis of the former is the pair- wise product of a basis of the latter). 

Suppose that K ^ const, so that at least one of the functions K^, Ky is not 
zero. We add the equation £'5 = and get a system £' C J^(2,3) of formal 
codimension 5. 

Its symbols gi C S^T* ® have dim^Q = 3, dimg[ = 1, dim 52 = and 
thus the only non-zero second (5-cohomology are iJ"'^(£') ~ M}, H^'^{£") ~ R^. 
There are two obstructions to compatibility - Weyl tensors W[ and W^2- The 
tensor W[ is proportional to 



Dividing this by bKy and simplifying modulo Ei, E2, E^, E4, E^^ we obtain the 
following expression: 

Eg = -35 \g\ hb ■ Vx + Qi ■ [u - w) + Q2 ■ v, 



Ux + \xU + XyV = 0, 

2vy -\- Wx + XxU + XyV = 0, 



Uy + 2Vx + XxV + XyW — 0, 
Wy + XxV + XyW — 0. 



-E5 = - [El ,E2,E3, E4]£ = 0. 
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where Qi, Q2 are certain differential expression of 5th order in the coefficients 
of the metric. The coefficients of -Eg, as well as others Ei, are not invariant, but 
the condition of their vanishing is invariant, and so can be expressed in terms 
of differential invariants. Indeed, 

Qi^J5- I^'V\9\{-K.){l3V\9\-Kl), 
Q2 = h- I^'^\{hM-Kl)"\hV\9\-'iKl), 

where 

^5 = (/sa - he) + [ha - hc){hc - &ha) - ^^hll 

is a differential invariant. Thus the equations Qi — Q2 = are equivalent to 
one condition J5 = 0. 

It is possible to show that this condition together with lib — implies 
hd = 0, which gives another proof of Darboux theorem [D] (proof in [Ko| ) that 
a Riemannian surface with 4 quadratic integrals is a surface of revolution. 

The second obstruction to existence of 4 integrals - tensor W2 - can be 
calculated similarly. Its vanishing is given by a scalar differential invariant of 
order 6 in metric, but it can be simplified modulo the conditions I^h = I^d = 
J5 = to the following expression: 

J4 = 3(/4a - hc)iha + ^hc)hc " l^hlliha + ^hc) + 25/|. 

Thus we obtain the following statement: 

Theorem 19. The condition of exactly 4 quadratic integrals can be expressed 
as 3 differential conditions on the metric: I^b = 0, J5 — 0, J4 — 0. 

If the compatibility condition EQ{mod Ei, E2, E3, E4, E^) = is satisfied, 
then the system £' is integrable. Otherwise we add this new equation and 
get (again in generic case, when the corresponding matrix of coefficients of 
derivatives is non-degenerate) the system £" with symbol 5" = 0, i.e. it is of 
Frobenius type. 

Its Spencer cohomology group H^'^{£") ~ R"^, so the obstruction to inte- 
grability - curvature tensor - has 3 components, represented by 3 linear 
equations relations on J*'(2, 3): 

E-jj = {u-w)+ A2jV = 0, j = 1, 2, 3. 

The expressions Aij are not invariant, but their vanishing is invariant and can 
be expressed via four differential invariants of order 6 of the metric; 

Jek = hk - P{l2,l3,hi,hj), k ^ a,b,c,d, 

where P is a quadratic function in I^j with rational coefficients in other variables 
(note that lee does not enter the formulae). All the expressions are rather long 
and shall be provided elsewhere. Let us indicate only equation Jga = 0: 
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ha = {TOOllhb - 825/2/|/5fc + 50/2/|/46(31/4a - IShc) 

+ 6hb{ha - hc)ieiL + 49/4^ - 37hahc + Qll) - 25/2/5b(-8/5a + he) 

- bhmlahb - 27hbll + 2hbhci-nha + 46/5c) 

+ ha{-i3hahb - 21hbhc + Shbhc) + 7/4^(4/56 - ll/5d))) ■ 

Theorem 20. The condition of exactly 3 quadratic integrals can he expressed 
as 4 differential conditions on the metric: ha = hb — he =^ hd ~ 0. 

Finally if E'jj are non-zero, we add these equations to the system. Com- 
patibility condition of the new system E'" are {Ejj)x = 0, {E'jj)y = 0, when 
expressed as linear functions on J°(2, 3) via the system £": 

Esi = Biiu + B21V + B31W = 0. 

Consider the matrix of coefficients of equations Eij,Eg,i: U = U{A,B). It 
always satisfies the condition rank(C/) < 3, because H is an integral of the 
geodesic flow. Also rank(?7) > if conditions of Theorem EOl are not fulfilled. 

Thus we have only two possibilities: If rankC/ = 2, the flow does not pos- 
sesses local quadratic in momenta integrals. Otherwise Esiimod Ejj) = 0. This 
means rank U = I and expressing this condition in differential invariants we 
get one condition of order 6 and four conditions of order 7 in coefficients of the 
metric (these long expressions will be omitted; note though that the above 4 
scalar differential invariants of order 7 involve 6 basic invariants hk, but the 
last of them hf does not enter the formulae): h = 0, hi = 0. 

These conditions give solvability of the system £, which yields us a 2- 
dimensional linear space of solutions generated by H and / - an independent 
integral of degree 2. We set □ — {h, ha, hb, he, hd)- 

Denote by S the singular locus of functions X{x,y), that are non-generic 
w.r.t. at least one one of the above steps (it consists of functions of one variable 
- metrics with non-zero Killing fields, and certain finite-dimensional families), 
which corresponds to zero denominators in hi- 

Let also be the set of A corresponding to the metrics with more than one 
additional quadratic integral (constant curvature or the conditions of Theorems 
[m and HOI)- We have proved: 

Theorem 21. There exists a polynomial vector-valued differential operator of 
order 7 □ : Cf^^iM.^) Cf^^{R'^;M.'^) and a residual subset S C C^^O^'^) such 
that for A G C^^{M'^) \ S the metric e^^^'^^ds'^^^^^ is Liouville (has quadratically 
integrable geodesic flow) iff X Cz £+ or □(A) = 0. Moreover, Liouville metrics 
in S are residual among all Liouville metrics. 

The singular locus S in the space of germs of Riemannian metrics (which we 
identified with C['^^(R^) only for convenience) is given by the condition lib = 0. 
The expressions of the above invariants together with a more detailed argumen- 
tation have appeared now in IK2I . 
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6.3. Gaussian curvature of minimal surfaces 

Consider a minimal surface Af ^ C M'' . The Gauss map defined on it depends 
on the curvature function and this function is unrestricted (i.e. can be arbitrary 
in a certain open domain; of course, it is non-positive, but there are no equality- 
restrictions) if the surface is considered abstractly (non-parametrized). But it's 
quite known that the Gauss map is not arbitrary, which is manifested by the 
fact, that the Gaussian curvature on the immersed (parametrized) surface is 
not arbitrary. We will describe precisely, which functions K on C M.^ are 
realized locally. 

u u —u^ 

So let be given as the graph z — u{x,y). Then Vi(m) — ^j'^^a'^^a^a 
is the Gaussian curvature operator and V2(u) = 

is the operator of mean curvature. Let = {u £ C^^{M.-^) \ V2('u) = 0} be 
the sheaf of minimal surfaces. We define Vi : H°° — > C^^{R'^) and denote the 
image by IC°°. Now we want to resolve this term: 

Theorem 22. There exists an algebraic differential operator □ : C^^{M.'^) — s- 
Ci'^p(IR^;M^) of order 4 md a finite- dimensional stratified suhmanifold S C 
Ci~ (R2) with K:°°\S = Ker(n). 

More precisely, there exist 4 polynomials ^Fj , _Fe , i^y , on the plane with 
coefficients depending differentially on K such that the function is realized as 
the curvature of a minimal surface iff they have a common root. 

The form of the operator □ = (Di, ^2) will be clear from the proof, though 
we suppress the formulas because of their size. Since the operator has singulari- 
ties, we obtain cases and K £ /C°° \5 iff Oi{K) = D2{K) = 0, while description 
of /C°° n 5 is given by some other operators, which we omit. 

Proof. Let us find the solvability criterion for the system Vi(u) = K{x, y), 
V2(w) = 0. Denote 

Fl = UxxUyy-uly-K-{l + ul + ul)'^,F2 = {l + ul)Ua:x-2UxUyUxy + {l + Ux)'^Uyy. 

The Mayer bracket of these operators is 

Fs = [Fl, F2] = aiiul + 2ai2UxUy + 022^^ + b, 

where an = K^i\n\K\)yy - ai2 = -K\\n\K\)xy, 022 = K^\n\K\)xx - 

AK^, b = K'^{{ln\K\)xa; + (In I^TD^j^) - 4K^. Notice that F3 has order 1, while 
generically the bracket of 2 second order operators after reduction is also of 2nd 
order. Therefore the system is compatible if F3 = 0, which is equivalent to 
= 0, i.e. the surface is a plane. 

Denote F4 = V^iFs), F5 = Vy{Fs). The system F2 = F4 = F5 = has the 
form: 

(l+uD -2UxUy (1+Wx) 

2(aiiUx + ai2Uy) 2(ai2Ux + a22Uy) 

2(aiiUx + ai2Uy) 2{ai2Ux + a22Wj/)y 









^xy 




bi 


_Uyy_ 




b2 
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where bi = -h^ - (aii)^u^ - 2{ai2)xUxUy - {a22)xul, 62 = "by ~ {aii)yul ~ 

2{ai2)yUxUy - {a22)yul. 

Resolving this for the second derivatives and substituting to Fi we get after 
multiplication by the square of the determinant of the above matrix a first 
order polynomial differential operator Fg. It has degree 8 by variables Ux,Uy, 
but reduction due to the system is of degree 6: 

Fe = -F^ + 4K{1 + ul+ ulf{2bFs - F|) 
= bl+bl+ib2Ux~biUy)^+4:K{l+ul+ul)^ {b^ +{aiiUx+ai2Uy)'^ + {ai2Ux+a22Uy)'^) . 

Denote — [E3,Eq] the Mayer bracket. It is the third 1st order PDE, 
which can be equivalently written as 

diF3,Fe) ^ diF3,Fe) 
d{x,Ux) d{y,Uy) 

Let us prolong the equation F^ = 0: VxiFe) = 0, 'Dy{FQ) = 0. Resolving 
these equations coupled with F4 — F5 — with respect to second derivatives 
(F7 = guarantees compatibility) and substituting to F2 = we get: 

2,d{F3,Fe) , ^, ^, fdiF3,Fe) diF3,Fe)\ ,^ ,^ ,25(^3,^^6) 



^ ^ d{x,Uy) '^\d{x,ux) d{y,uy) J ' d{y,Ux) 

Thus we obtain four 1st order PDEs F3 = 0, Fq = 0, F7 = 0, Fg = 0, which are 
polynomials by Ux,Uy of degrees 2, 6, 7, 9 (they do not depend on u) and are 
differential operators by K of orders 2, 3, 4, 4 respectively. 

Solvability of the system Fi = F2 = is equivalent to the claim that poly- 
nomial by Ux,Uy system F^ ^ Fq = Fj = Fg = has a solution. The latter is 
equivalent to 2 conditions Di = 0, ^2 = 0, algebraic by the corresponding coef- 
ficients, which are differential operators of K. Therefore we set □ = (Di, 02). 

The set S is formed by functions K for which some of ^3,^6,^7,^8 be- 
come dependent or singular. This set is given by a collection of overdetermined 
systems of PDEs of finite type and hence is stratified finite-dimensional. 

Let us demonstrate a geometric idea behind this proof. Fi and F2 are Monge- 
Ampere operators on the plane, which means ( IL2I ) that they are given by 2- 
forms f2i,J72 on J^(M^): Equations Fi{u) = can be rewritten as = 0, 

where ji{u) C J^(M^) is the jet-section determined by u and 

fli — dux A duy — K[\ + + Uy)^dx A dy, 
fi2 = (1 + Uy)dux /\ dy + UxUy[dx A dux + duy A dy) + (1 + u^)dx A duy. 

Since fi^ do not depend on m, the construction descends onto r*R^ — J^(M^)/M^, 
where we have in addition the canonical symplectic 2-form fip. 

We search for a Lagrangian surface, which is isotropic w.r.t. ili, 0,2 and lies 
in the hypersurface — {F-^ = 0} C r*R^. Let £,i be the kernels of the 2- forms 
f^i restricted to S^. They fail to be in general position only along a surface 

= {F^ = Fq = 0}, which is equivalently determined by the equation 

E2 ^ {a; e I rank(^o(a;),ei(a;),e2(x)) < 3}. 
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This surface is ramified over with the fiber consisting of no more than 12 
points (precisely this number if counted over C and with multiplicities). 

Local solvability means that over a neighborhood in M? at least one compo- 
nent of is isotropic for all fli. By characteristic property of this surface 
we need to require only two restrictions riils^ to vanish. For instance, we take 
the conditions ^0,^2 C TE^ . They are given by the equations QqA dF^ A cIFq = 
and fl2 A dFs A dF^ = 0, which correspond to the operators F7 and Fg. □ 

Example 3. Consider the family K = A;e"^+^^. Then it realizes the curvature 

of a minimal surface iff k = 0. In fact, the quadric F3 — b is definite and 
diib > 0. So F^k < and the inequality is strict if k ^ 0. 

Notice that ^3 is non-singular UK ^ 0. The set /C°° is a subset of the 
solutions to F3 = Fq = 0. Let us describe non-holonomic solutions to this 
system (ignoring the compatibility condition F7 = 0). 

Let's parameterize the quadric F3 = 0: Consider a line = t,Uy = \t. It 
meets F3 = at the points with t"^ = —b/p{X), p{X) = an + 2ai2A + a22A^. 
Substituting this into p{X)^Fq = we get a polynomial of degree 6: 

+ 4K [p{X) - 6(1 + A^)] ' [6V(A) - 6((aii + aisA)^ + (ai2 + a22A)2)] . 
Thus {F3 = 0} n {Fg = 0} corresponds (2-to-l) to the roots of qe{X) = 0. 

Example 4. Consider the family K = ip{x). For an open set of such functions 

the equation q(,{X) = has 6 positive roots and so the system F^ = Fq = 
has 12 non-holonomic solutions Ux,Uy. None of them satisfies the other two 
equations F^ = Fs = 0, save for the case Ux = consti, Uy = const2. Actually, 
the solutions depend on x only and so Uxy = Uyy = 0, whence Uxx = and 
K ^0. 

Lot So C S he the set of functions K, such that {F3 = 0} = {Fq =0}. It is 
given by the condition g6(A) = 0, which is a system of 7 third order PDEs on 
K. As a by-product of the proof we get the following statement: 

Proposition 23. For K E /C°° \ Sq the set of minimal surfaces through the 
origin with this curvature Vi^{K) n 7i°° n {m(0, 0) — 0} has cardinality at most 
12, though generically this number is 2 (respectively 6 and 1 if we the function 
u is considered up to the sign). 

Proof. In fact, all the equations are symmetric to the change u 1— > —u. 

The sub-system F^ = Ff, = Q with respect to Ux, Uy has as maximum ord(F3) • 
ord(F6) = 12 algebraic solutions. Some of them may not satisfy the other 
constraints Fj = Fs = and generically (in /C°° \ <So) only 2 do satisfy. □ 

The above statement does not hold for K = G Sq. But we suggest there 
are no other examples. To see the reason let us consider the overdetermined 
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system of 7 polynomial differential equations on K of order 3: (?6(0) = 0, (7g(0) = 
0, . . . ,96^''(0) = 0. One can expect that it has no other solutions than K = Q, 
but this is not so. 

For instance, qe{\) = follows from b = Q and "^^^ = 0. This latter is 

0-12 ^22 

equivalent to the system detHess(ln \K\) =0, A In \K\ = 4i^, which though in- 
compatible is solvable and has the solution K = (p{xzLy), where ip' = 2(py^Lp + c. 

But for this choice of K the other two equations = Fg = are not 
satisfied unless (p = 0. In fact, then F^ — and Fq — are both equivalent to 
Ux + Uy — 0, which coupled with Fi = F2 = gives K = 0. 

The computer programs do not give us other solutions to the above system 
of 7 third order PDEs on K and so we conjecture that 

/C°° n 5o = 0. 

Remark 15. This our conjecture that except the plane case a minimal surface 
is restored (up to ±) from its Gauss image in maximal 6 different ways is sim- 
ilar to the known Gronwall conjecture about webs on the plane. It says that a 
linearizable 3-web has maximally 11 linearizations (and generically only 1), cf. 
\GLi^ . Both problems are basically algebraic and have equal complexities. 

Similarly one investigates the problem, when two functions K and H on a 
parametrized surface can be realized as Gaussian and mean curvatures, with 
the surface realized as a graph (projection R'^ yields the parametrization) . 

This is an analog of the classical Bonnet problem of realizing two quadrics 
as the first and the second quadratic forms on a surface. Bonnet theorem states 
that compatibility and solvability of this problem is the system of one Gauss 
and two Kodazzi equations. 

For realization of the curvatures K, H compatibility is equivalent to the 
condition K = H = 0, i.e. the surface is plane (notice that the solutions space 
is 3-dimensional, not 4-dimensional as one can expect after ^ 35.41 but this is due 
to non-genericity of the condition). Solvability leads to an operator of order 4, 
similar to the above □. Thus we get solution to generalized Bonnet problem. 



6.4. Quantum integration 

Consider the algebra of scalar linear differential operators A = Diff (1, 1) on 
the manifold M filtered by the C°°(M)-modules Ak of order < k differential 
operators. Let V = (Bk>o'Pk = g^{A) be the corresponding graded module. Here 
Vk = Ak/Ak-i consists of degree k homogeneous in momenta polynomials on 
T*M. Thus V = SV = ®S''V, where V is the C°°(M)-module of vector fields 
on M. The canonical Poisson structure on V is given by the formula 

{VimodA, V2mody4,;} = [Vi, V2] mod A+z-i, 

where the bracket in the r.h.s. is the usual commutator (or Jacobi bracket). 
In other words, the mapping a = smbl : A ^ V is a homomorphism of Lie 
algebras. 
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By quantization one understands an inverse map q^, -.V ^ A, i.e. a collection 
of morphisms qk '■ Vk Ak splitting the sequence 

^ Difrfe_i(l, 1) ^ Difffc(l, 1) ^ S'^V ^ 0. 

This map allows to introduce a new non-commutative associative product 
on V: 

a-kb^ g7^(g,(a) o q^{b)). 

These kinds of products are important in the deformation quantization. Moyal 
|Mo| and other star-products IL3I are obtained by specifying the morphism q,. 

Denoting by pk the homogeneous Pfc-component of a polynomial p V we 
observe the relation 

{a, 5} ^ {a-kb- bka)k+i-i = smh\k+i-ii[qk{a), qi{b)]), a <E Vk,b e Vi 

between the commutator, Poisson bracket and the star-product. 

Consider a mechanical system with a Hamiltonian h € V. Its quantization is 
given hy H = If the choice of is subject to certain connections (IL3I), 

then Riemannian metric produces the Laplace operator, choice of potential - 
Shrodinger operator etc. 

Suppose the classical system is integrable in Liouville sense, i.e. there exist 
functions /i = h, f2, ■ ■ ■ , fn & V {n = dimM) functionally independent a.e. 
which Poisson-commute {/i, fj} = 0. We wish to quantize this picture. 

Definition 6. Differential operator H is called quantum completely integrable 
if there exist commuting differential operators Fi = iJ, , . . . , G A, which 
are independent a.e. 

Clearly then the system h = smbl(_H') is Liouville-integrable with integrals 
fi = smbl(_Fi). The quantization poses the inverse problem: To find quantum 
integrable system {H, Fi) by the given classical (ft, fi). 

This problem was solved for many classically integrable Hamiltonian systems 
( |Pe| ) using different approaches: analytical, Dunkl's differential-difference op- 
erator [Duj . Moyal quantization, via geodesic equivalence |MT] and others. 

We discuss one of them, which is closely related to our integration method. 
It was proposed in |He| and is based on universal enveloping algebras. 

Consider the rigid body equations, which is the Hamiltonian system on 
T*SO{3) with Hamiltonian h = \ Y.T~^Pi + E 

7iXi, where the base 

coordinates and pi are the corresponding momenta. 

Let Xi = q*{xi) at i? € S0(3) be equal to {Rci, e), where is an orthonormal 
basis and e some unit vector in K^, and Pi = q*{pi) be the left-invariant fields 
exp{Ei) generated by the basis Ei S so(3) given by the relations Ei{ei) = 0, 
Ei{ei±i) = iciipi, i G Z3. Then the subalgebra of ^(S0(3)) generated by Pi, Xi 
is isomorphic to the universal enveloping algebra of so(3) k M'^: [P^, Pj] — eijkPk, 
[Pi,Xj] = djkXk, [Xi,Xj] = 0. 

The quantized Hamiltonian has the form H = q^{h) = '''i'^^? + J2 liXi 
and we denote it also by Fi. We have two Casimir functions in J7(so(3) k M.^): 
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F2 = ^3 = J^-^i-^i- To achieve complete quantization of all known 

integrable Euler equations we must quantize the forth integral. It is as follows: 

Euler case: 7, = 0. Then F4 = i^^^ 

Lagrange case: ti = r2, 71 = 72 = 0. Then F4 = P3. 

Kovalevskaya case: ti ^ T2 ~ T3/2, 73 = 0. Then F4 — KK + KK — (74, 
where K ^ ti{Pi + iPj)^ - 2(71 + ij2)iXi + 1X2) and .94 = StI{P^ + In 
all these cases [H, F4] = 0. 

Goryachev-Chaplygin case (conditional integrability): ti — T2 = 
73 - 0. Then F4 - ri(p2 + p2)p^ _ X^ijiPi + 72^2) - 53, where 
93 = \il2Xi - J1X2) + JT1P3. In this case [H, F4] = Ti{j2Pi - 71-^2)^3- 

Note that according to ||L3j any quantization is determined by two linear 
connections (electromagnetic field and gravity) and a collection of tensors gk ■ 
S^V — !■ Ak-2- The first two of the above integrable cases are obtained from 
the classical scheme via the operator and trivial tensors 52 and gi. In the 
Kovalevskaya case one needs a second order correction 174 and in the Goryachev- 
Chaplygin case a first order g^ (these corrections were found previously in jHe] , 
but the explanation of actual orders meaning was lacking). 

Theorem 24. In each of the classical integrable cases the obtained quantum 
integrals allow to integrate the Shrddinger operator L[u\ — Ut — H{u) — AiU 
(where Ai is the spectral parameter) classically: The system 

L[u\ = 0, F2{u) = X2U, F3{u) — A3U, F4{u) = A4M 

is of finite type and compatible (A3 — in the Goryachev-Chaplygin case). 

Proof. In the first three cases, where we have a commutative collection of 
integrals Fi, F2, F3, F4, the statement is rather known. However, the Goryachev- 
Chaplygin case seems to be quantization by analogy and its meaning is given 
by the above statement, which follows from our compatibility criterion. □ 

Remark 16. The above result is constructive, not mere an existence statement. 
In fact, the considered system is of Frobenius type and one reduces its integration 
to a certain system of ODEs, which may be integrated via symmetry approach. 

So far we obtained only local solutions, but if we are given a global Shrddinger 
equation on a manifold M , then we get a topological restriction for a compat- 
ibility/solvability. Namely the monodromy operator determines which spectral 
parameters A are admissible. For them and only for them we get a closed leaf 
of the foliation corresponding to the above Frobenius system. 

Note that in the definition of quantum integrability we used smbl-map, not 
. In fact, q^^{Fi) need not to commute, only their top-components. Denote 
bypfe -V the natural projection. Then commutation of gfc (a), qi{h) implies 

that {a, 6} = pk+i-iq'^^ {[q*{a) , qi,{b)]) = 0. So the quantization philosophy 
suggests to make the construction so that the functions commute w.r.t. the 
deformed bracket 

{a,b}q ^ q:^^{[q*{a),q^{b)]). 
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This is equivalent to the quantum integrability problem: given a system 
/i,...,/„ involutive w.r.t. the Poisson bracket describe the quantizations 
such that the set is still involutive w.r.t. the new bracket {, 

However our compatibility result suggest that it is equally important to 
search for deformations giving sub-algebras, i.e. {fi,fj}q = ^c'^jfk- In the 
classical case they correspond to invariant submanifolds of the Hamiltonian sys- 
tem. Thus our main result interprets as quantization of conditional integrability 
in the classical mechanics. 
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